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SYNOPSIS 

A mathematical description of the process units is 
necessary for the design, optimization and control of a 
chemical plant. These mathematical relationships are often 
obtained by regression analysis of historical plant data. 

The reliability of these regression models is limited by 
system noise. Also these models are valid only for the range 
of the variables in which data were collected. 

A better model can be obtained by considering all the 
physical and chemical processes taking place in any differentia 
section of the unit under consideration. This work describes 
process simulation of ammonia plant wherein, for each important 
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urLii; sucli nis.‘tii6iii9,t;jLca-l mod-els liave been developed. 

Ammonia, wMcb. is a basic fertilizer ingredient, is 
industrially synthesized from nitrogen and hydrogen. The 
process consists of steam reformation of a hydrocarbon 
(natural gas or naphthas), conversion of CO in the reformed 
product to COg, removal of CO 2 hy absorption, conversion of 
any residual carbon oxides to methane and finally synthesis 
of ammonia from this gas which contains '^ 2 * ^ 2 ^ 

The important process units are primary and secondary reforme] 
high and low temperatiire shift reactors and ammonia synthesis 
reactor. 

In the primarjr reformer model development, a method to 
calculate transfer of heat to reformer tubes from flames in 
addition to those from flue gases in the fire box, has been 
developed. The total transfer of heat, from flames and flue 
gases, to a differential section of a reformer tube is 
considered simultaneous with the transfer of heat to the 
reacting gases and the reaction taking place inside it. Erneti 
rate expressions describii]g steam -hydro carbon reactions, under 
varying conditions of pressure, have been obtained by modifi- 
cations in the available ones. This model has been used to 
check performance of reformers of different sizes (reformer 
tube heated length -3.28 to 12.2 meters), using different 
hydrocarbon feeds (natural gas or naphthas) and operating 
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•under different pressures (13.7 to 35.0Kg/cin3 and feed 
temperatures (371 to 456 , 4 '^G). Agreement between plant data 
and sim'ulation results is generally very good. 

In the secondary reformer model, the temperatxare and 
composition at the inlot, duo to mixing of air with high 
temperature product from primary reformer, is calculated, Rato 
of refoimation reaction is calculated by considering diffusion 
of methane to the catalyst pellet surface to bo the limiting 
factor. Match of the calculated inlet temperature, amount of 
air used and the outlet tempera buro and composition i\rith the 
respective plant data confirms the validity of the method of 
calculation. 

Rate equation for -water-gas shift reaction over a high 
temperature shift catalyst (Pe 20 ^-Cr 20 ^) has been developed. 
This rate equation takes into account the effects of tempera- 
ture, pressure, ago of the catalyst, H 2 S content of tho 
reacting gases as well as the intraparticle diffusional resis- 
tance in the catalyst pores. It has been used in the high- 
temporature sliLft reactor model to check the performance of 
reactors having single or multiple beds and working undor 
different conditions of temperature (603-744^K) , pross-ure 
(16,4 to 30 Kg/cm^), composition of the reacting system (CO - 
9.0 to 40,6 per cent), catalyst age (60 to 560 days), H 2 S 
concentration (less than 0,5 to 60 ppm) and steam to carbon 
monoxide ratio (5 to 21.7), Tho close agreement between the 



plant data and predicted values validated the model, 

For lo'w-tempcraturo sliift catalyst (CuO-ZnO), the 
effects of temperature, pressure, ago of the catalyst and 
intraparticle diffusional resistance on the rate of reaction 
arc accounted along the same lines as for high temperature 
catalyst. Performance of reactors using catalysts of differor 
ago (60 to 560 days) and working under different conditions 
of pressure (l6.0 to 28.9 Kg/cm^) and temperature (183-228*^0) 
can satisfactorily be predicted with this model. However, 
the change in composition as well as temperature is rather lo'u 
in this reactor and hence the significance of the model is 
not obvious. 

A workable method to calculate the effectiveness factor, 
for a complex reaction, by partially solving the transport 
equations has been developed. Appropriate rate equation 
describing synthesis of ammonia over an iron catalyst has 
been selected. This, along with the method developed to 
calculate effectiveness factor, has been used in a mathomatica 
model to describe the ammonia synthesis reactor performance. 
Both adiabatic as well as non-adiabatic beds have been consi- 
dered, The model very successfully ropresenbs the plant 
performance. 

These unit models arc finally unified in an overall 
plant simulation which assumes ideal performance for the 
units which have not been separately considered. This program 



is used to ch-cck tho porf ormanc o of two existing ainmonia 
plants each of 600 tons/day capacity, working with, natural 
gas and naphtha as tho food stock and at difforent conditions 
of oporation. The agreement between plant performance and 
simulation results have been chocked at various points 
throughout the plant. 

The simulation models have been used only to chock 
the performance of existing units or plants* However, in 
view of their validity over diverse conditions of oporation, 
those models could as well be used for design purposes. 



CHAPTER 1 


IHTRODUCTIOH 

With increase in world fertilizer demand the need 
for the basic fertilizer ingradient, ammonia, has also 
increased. The design development in mid sixties, of 
single train ammonia plants employing centrifugal compressors, 
facilitated large scale production capacity. As a result, 
numerous plants have been constructed around the world which 
produce 600-1500 tons of ammonia per day. The plants use 
feed stock ranging from natural gas to naphthas in the 
boiling range upto 220°C. 

Figure 1.1 shows a simplified process flow diagram 
for such an ammonia plant. After being mixed with steam, 
the desulphurised process hydrocarbon enters a Primary Reformer 
(PR), The outlet stream containing 00, CO 2 » and steam, 
is mixed with air at the Secondary Reformer (SR) inlet. The 
quantity of air used is such that after complete conversion 
of CH^ and CO to CO^ the hydrogen to nitrogen ratio 
in the process gas is slightly more than three. The high 
temperature SR product is cooled (Heat recovery) and sent 
to High Temperature shift reactor (HT) where GO reacts with 
steam to yield additional CO 2 and H^. In the Low Temperature 
shift reactor (LT) that follows, the GO concentration in the 
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gas is reduced to a very low value closely approaching 
equilibrium at a lower temperature. Carbon dioxide is 
subsequently removed from the process gas in 00^ absorber 
and remaining carbon oxides are converted to methane in the 
methanator. The methanator exit i.e. synthesis gas, is 
compressed and mixed with the recycle gas to make feed for 
the ammonia synthesis reactor. Reaction in this reactor, 
between nitrogen and hydrogen of the synthesis gas, yields 
ammonia which is the final product. Ammonia is separated 
from the product stream, a part of the remaining gas is 
purged to maintain low concentration of inerts in synthesis 
gas and the rest of it is recycled. 

The general engineering philosophy behind the single 
tram ammonia plants is to employ higher pressures in the 
reforming section (PR and SR) and lower pressure in ammonia 
synthesis section, compared to those in older plants which 
use reciprocating compressors. The centrifugal compressors 
used in these plants are dnven by steam turbines. Required 
steam is generated by the utilization of sensible heat from 
the product gases of the synthesis reactor and the SR as 
well as waste heat from the PR flue gas. While a major 
advantage of this design is to obtain an extremely efficient 
utilization of the heat generated in the plant, its operation 
is highly sensitive to both process and equipment performance. 
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This ■work takes only the process performance into consideration 

A mathematical description of all important process 
units taking into consideration all the important physical 
and chemical processes taking place inside them is necessary 
for a successful design and operation of the plant. More- 
over, when a computer is used for control of a chemical 
process, a mathematical description of the process units is 
necessary to perform such tasks as feed for-ward and feed 

hack control, as well as off line and on line optimization. 

1 

Historical plant data are often used to obtain such a mathe- 
matical description by means of regression analysis. However, 
reliability of these regression models is limited by system 
noise. Also these models are valid only for the range in 
which the values of the variables were obtained and extra- 
polation IS generally not very useful. 

Availability of a priori knowledge of the behaviour 
of the critical units greatly facilitates the design of 
computer control strategy of a plant. Off line simulation 
of the actual plant units by means of analog or digital 
computer can provide such a knowledge, provided that the 
mathematical models chosen for simulation describes the plant 
with reasonable accuracy. Once the accuracy and generality 
of the plant simulation in predicting the plant behavio'ur is 
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established, the same models can be used for designing new 
plants. Also this could be reduced by various simplification 
techniques for use in on-line control and optimization. 

The important process units in the ammonia plant are 
the primary and secondary reformers, the high and low 
temperatTire shift reactors and the ammonia synthesis reactor. 
Simulation model for each of these units is described in 
this work. These unit models have been tested under diverse 
conditions of operation and finally a unified computer program 
using all these models has been prepared to check the per- 
formance of two ammonia plants. Each of these are of 600 tons/ 
day capacity but one of the plants is based on natural gas 
while the other uses straight run naphtha in the boiling 
range 200 to 210^0, as feed stock. ' 

Chapter 2 describes the models for primary and secondary 
reformers. In Chapter 3 model development of high and low 
temperature shift reactors is discussed. Model for ammonia 
synthesis reactor is discussed in Chapter 4, The 5th Chapter 
deals with t unification of all these unit models in an overall 
simulation program for ammonia plant. 
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) % 

SIMULATION OF S^^EAM-HTDROCARBOR REFOUH^S 

Steam reformation of hydrocarbons is carried out 
in two consecutive reactors, namely, Primary Reformer (PR) 
and Secondary Reformer (SR). PR consists of rectangular 
furnaces in parallel. In each furnace, catalyst tubes 
are arranged in two rows parallel to the two walls in 
which a large number of burners are embedded (Figure 2.1), 
Reactions over the latest steam -reforming catalysts are 
so fast that thermodynamic equilibrium, corresponding to the 
process temperature, is closely approached. Since higher 
temperature favours low methane content in the process gas, 
it is desirable to attain high PR exit temperature. However, 
limitations of heat transfer at higher temperature through 
metallic tubes put a limit to which the temperature can be 
raised. 

Process condition close to equilibrium at higher 
temperature is attained in SR. Air is added to the SR inlet 
stream (PR exit) to provide nitrogen for ammonia production 
and the necessary heat, obtained by the combustion, for 
conversion of the res3 "'ual methane in PR exit. There is no 
heat transfer through the reactor (SR) walls and performance 
IS dependent on adiabatic partial combustion of inlet stream 
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followed "by reforming over the SR (nickel) catalyst. 


2.1 PRIMARY REFORMER 

The feed to PR (hydrocarbon-steam mixture) enters 
the catalyst packed reformer tubes at the top. Ehdothermic 
heat of reaction as well as heat to increase process gas 
temperature is supplied by large number of burners distri- 
buted uniformly over the walls (Figure 2.1). 

There are some works reported in literature which 
deal with the problem of heat transfer in top-fired [27,30] 
and side fired [16] steam -hydrocarbon reformers. Some 
others have simulated steam hydrocarbon reactions in a 
reformer [23,24]. The studies related to heat transfer in 
a side-fired reformer consider only the total amount of 
heat transferred over the entire length of reformer tubes. 

In the studies on simulation of reformer reactions, tube 
wall temperature profile has been assumed. The objective here 
IS to find the conditions of the process stream at any point 
inside the reformer tube given the inlet conditions of 
process feed and burners. 

For this purpose, kinetic rate expressions descri- 
bing reactions inside a reformer tube have been developed. 
Transfer of heat from flames and flue gas to reformer tubes 
has also been modelled. With the help of these, a mathematical 



model, taking into account the processes inside a differen- 
tial reformer tube section and transfer of heat to it by 
flames and flue gas, is developed. Calculation program 
used for its integration provides the composition, temperature 
and pressure at each point inside the tubes and at the exit. 

It also gives the temperatures at all points on outside and 
inside surfaces of the reformer tubes and the ^ount of heat 
transferred through these surfaces to the process gas. 

2.1.1 RATE EQUATIONS 

The steam -reforming process generally uses nickel 
catalysts and feedstocks ranging from natural gas to naphthas 
of final boiling point 220° 0. The chemical reactions that 
take place in catalytic reforming of these hydrocarbons are 
exceedingly numerous. Particularly in case of naphthas, 

•which have very large number of components, it is practically 
impossible to determine the mechanisms and kinetic rate 
expressions of different reactions taking place. Therefore, 
many speculations have been made regarding the mechanism of 
steam -naphtha reaction and consequently to the development 
of models for the same. The model generally proposed, used 
in this work and elsewhere [23,24] is based on refoimation 
experiments carried out with hydrocarbons heavier than methane 
e.g. propane, butane, hexane and benzene [6,12,32]. Prom 
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such, experiments it is deduced that methane is the primary- 
product of reaction "between steam and heavier hydrocarbons. 
Accordingly, in the present model, methanation is taken to 
be completed at the entrance to the reformer i.e. all 
hydrocarbons heavier than methane are assumed to hydrocrack 
to methane. Consequently, the reaction system inside a 
refoirmer tube is described by the kinetic rate expressions 
for steam-methane reaction, irrespective of the process feed 
hydrocarbon. 

Any pair of the following four reversible reactions 
will account for the stoichiometry in steam -methane reforming 


CH4 + H^O CO + 3^2 (2.1) 

CO + H^O CO 2+ H2 (2.2) 

CH^ + 2H2O CO 2-^ 4H2 (2.5) 

OO2 + GH^ 2C0 + 2H2 (2.4) 


While the choice of equations is important when 
dealing with kinetic relationships, any pair of these given 
here may be selected when dealing with equilibrium relationship; 

Taking different pairs from among the four reversible 
reactions described above, maiy kinetic models for steam- 
methane reactions have been proposed [1,2,7,12,29]. These 
rate equations differ considerably from one another but 
there are some common features as well. Firstly, the rate of 
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disappearance of methane is first order with respect to its 
concentration for almost all the cases referred to, although 
the pair of reactions used to describe the process are 
different. Secondly the forward rate constant, is dependent 
on pressure. 

The differences among the equations, for rate of 
disappearance of methane, could be due to more than one 
reason. However, the most probable one seems to be the 
neglect of pore diffusional resistance, which is significant 
in large size catalyst pellets when the reaction is very 
fast. The inter-particle mass transfer resistance is insigni- 
ficant under the high mass velocity conditions existing in 
primary reformer. Moreover, in case of shift reaction 
(reaction 2, 2) , seemingly different rate equations obtained by 
various workers [3,9,10,18,22,25,28] could be reduced to a 
single (intrinsic) rate equation by accounting for pore 
diffusional resistance [31, 33], Based on these observations 
it has been assumed that pore diffusional resistance for 
reformation reactions is very high. Pair of reactions (2,l) 
and (2,2) has been taken to describe the reaction and the 
rate is evaluated by the following first order kinetic 
expressions [19] 

‘ Poo 

r^ = a^ exp(-Ea^/RT) (p^j^^- 


) (2.5) 
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1*2 — 3.2 * 6 xp (—232/^^^) ]^ * T) ^ (2*6) 

eq2 ^ 2 *^ 

where r^ and r 2 are rates of reactions ( 2 .l) and ( 2 . 2 ) 
respectively; a^ and are the pre-exponential factors, 
constant at given pressrire; Ea^ and Ea 2 are apparent energies 
of activation, and are equilibrium constants, P is the 
total pressure and p is the partial pressure pressure of 
component j . 

The values of apparent energies of activation and 
pre-exponential factors used in this work are 

Ea^ = 8,780 K. cal/kg.mole 

a^ = 127 Eg. mole/hr, Eg. of Catalyst, atm. 

Ea 2 = 13880 E. cal/Eg.mole 
a^ = exp(8.23) Eg. mole /hr. Eg. of catalyst, atm. 

The values of a^ and a 2 are at a total pressure of 1 
atmosphere. Values of Ba^ and a^ are due to Akers and Camp[l], 
Those of Ea 2 and a 2 have been taken to be the same as those 
over an iron oxide catalyst (shift reaction catalyst). Values 

of Ea 2 and 32 are from Singh and Saraf [33]. 

If the rate of a first order reversible reaction, 
with significant pore diffusion resistance, is expressed in 
terms of partial pressures of reacting components, the rate 
constant is inversely proportional to square root of the total 



13 


pressure [31]. So the values of pre -exponential factors, a^ 
and a^, at any pressure other than one atmosphere could be 
obtained by dividing its value at one atmosphere by square 
root of total pressure i.e. VP. Assuming the ideal gas 
conditions to hold and expressing the partial pressures as 
multiples of total pressure and mole fractions we have from 
equations (2.5) and (2.6) , 


127VP exp (-§^1^) (x 




‘CO. 


‘CH, 


K 


eq.1 




^2^ 


‘CO, 


') (2.7) 


r, = exp (- 


RT 


+ 8.23)V'P (xqq- £ 


eq2 H^O 


) ( 2 . 8 ) 


The above equations are valid under all operating 
conditions. 


2.1,2 HEAT TRANSPER 

The arrangement of tubes and burners alongwith the 
path of flue gas flow in a side fired reformer is shown in 
Pigure 2.1. The heat transfer problem in such reformers can 
be conveniently split into two parts, namely, the transfer 
of heat from flames and flue gas to the reformer tubes and 
from inside walls of the reformer tubes to the reacting gas 
mixture. 

Radiation is the prevailing mode for transfer of 
heat from flue gas and flames. Heat transferred by conduction 
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and convection is negligible compared to the total amount 
of heat received by the reformer tubes. Flue gas contributes 
larger part of heat in total transfer where as flames 
emanating from large ntmber of burners provide the rest. 

For the sake of ease in calculation, the flames and flue gas 
are assumed as two distinct sources. Radiative transfer 
from these two sources to the tubes is assumed jhot- -to interfer^ 
with one another and the sum of the two is the total heat 
transferred. Both these sources consist of a mixture of 

i 

radiating gases and the their emissivities enter all calculatio' 

! 

I 

Emissivities of Flames and Flue G-as ; | 

The flames and fully burnt combustion products (flue 
gas) of most conventional gas and oil fuels are real gases, 
since the radiation from them is mainly due to carbon dioxide 
and water vapour which emit in discrete bands. However, the 
form of the relations for radiative heat transfer in an 
enclosure is greatly simplified if the gray gas assumption is 
made. Therefore, it is advantageous that the mathematical 
formulation which characterises a gray-gas system is retained 
in the representation of a real gas emission. 

This IS achieved by representing the emissivity of 
a real gas as the sum of the weighted emissivities of a 
number of gray gases [20]. The emissivity of a gray gas is 
given by 
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Eg = 1 - exp (-k^ pL^) 

where k is the absorption coefficient, p is the partial 

O 

pressure of gas and is the mean beam length. The total 
emissivity of a real gas, in terms of weighted emissivities 
of a number of gray gases, is given by 


1,' ^g,n ( 2 - 9 ) 


where a and k are weighting and absorption coefficients 

for component n respectively and N is the total number of 

components. If the number of teims in the above summation 

i.e. N, IS very large, a may be considered as the fraction 

6 » ^ 

of energy in a blackbody spectral region in which the effective 

absorption coefficient is k . However, k and a may 

be more correctly considered as numbers which make the 

series in equation (2.9) fit the given function of e vs pL . 

g m 

Although in theory the emissivity of a gas approaches one as 
the value of pl^^^ becomes very large, at all practically 
attainable values of it the emissivity is considerably less 


than one. Therefore, in fitting the emissivities of real 

gases with finite number of terms, it is common to have one 

term with k =0. This physically corresponds to the windows 
s y ^ 

in the spectrum between strong emission bands and is sometimes 
referred to as the clear gas component. 
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In such a representation k is held constant and 

gy n 

tremperature dependence of emissivity is transferred to 
weighting coefficients by expressing a „ as 

Sy n 


a = + 

g,n l,n 


b^ T 
2,n g 


( 2 . 10 ) 


where b^ and b„ ^ are constants for each value of n and 
1 » n » n 

T is the gas temperature. In the above representation Taylor 

o 

and Foster [35] have used a 1 clear + 3 gray gas components 
model to fit the emissivities of combustion products of 
commonly employed gas and fuel oils as 


e ?= 

g 


4 

E 

n=l 


(b. 




T )' (1-e 


— k (p +p )l \ 


( 2 . 11 ) 


where p^ and p^ are partial pressures of carbon dioxide and 
water vapour respectively, in the radiating gas mixture 
(flue gas or flames). This correlation, with values for 
constants in the reference, have been used here. 


( 1 ) Transfer of Heat to Reformer Tubes ; 

( Q-) Ra-diative Transfer from Flue G-as ; 

The enclosing refractory walls of a reforming 
furnace (fire-box) do not transfer any heat except that lost 
to atmosphere. This loss is generally assumed to be equal to 
the amount of heat transferred to it by convection and 
conduction from flue gas. So, there is no net radiative 
exchange between the enclosing walls and radiating gas and 
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tubes. Such, a surface is termed a no flux surface. 

A uniform flue gas temperature is assumed. The 
burning gas emanates at considerable momentum, from the 
burners, leading to a thorough mixing of gases inside the 
furnace. Measurement of temperature at different points 
throughout a reforming furnace also showed an insignificant 
variation, justifying the assumption. 

It IS further assumed that reformer tubes are gray 
and the disposition of surfaces is such that the view factor 
of the tube (sink) surface from a point on the refractory 
walls IS same as that from any other point. The latter 
assumption implies that sink and walls are intimately mixed 
i.e. the walls are speckled. 

In view of all the assumptions made, "the problem can 
be visualized as that of radiative heat exchange between an 
isothermal gray gas and a single sink zone (tube surface), 
the enclosure being completed by refractory s urfaces grouped 
together into a single no flux zone. Prom the available 
derivation for such a case [20] the net exchange between flue 
gas and reformer tubes can be written as 


Qg.^t = ^®g“®t^ 

where 


1 



1 ) + 


1 

black 


( 2 . 12 ) 


(as^)j^ 


( 2 . 13 ) 
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A. 


black - 1 + egAl-eg)F|it ^ 


(2.14) 


E. 


S 


A A. 

crTg , E^ = ^ St ef an-Boltzmam 


constant, T and T aro the gas and outer tube surface 
g o 

temperatures respectively, is the total exchange area 

between flue gas and reformer tubes and black 

total exchange area when tube surfaces are assumed to be 
black, e and e. are emissivities of gas and tubes respectively, 
A^ and A-^ are the surface areas of tubes and refractory walls, 
and PRt IS the view factor of tubes with respect to refractory 
walls. can be determined from charted values [20, 26] 

by considering an imaginary surface parallel to the tubes 
(and just outside them) and multiplying the view factor between 
the refractory and this surface by that between the surface 
and tubes. 

It is clear from equations (2.13) and (2,14) that the 
total exchange area, (G’S^)j^» is independent of the tube 
surface temperature. In addition, the assumed intimate mixing * 
of reformer tube surface and that of the refractory wall makes 
it clear that there is no distinction between different 
elements of reformer tube surface due to geometrical disposition. 
Hence, the exchange area for a unit reformer tube surface is 
constant irrespective of its location and it is equal to the 
total exchange area divided by the total tube surface area. 

So the net exchange between furnace gas and an element of 
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unit area on a reformer tube can be written as (Equation (2.12) 


1 


g 







(2.15) 


where is the total number of reformer tubes, 

is tbe total outer tube surface of a reformer tube and D„ 

o 

and L are outside diameter and length of a reformer tube 
respectively. 


( b) Radiative Transfer from Flames ; 

Following _assum_ptions are made: 

1. Burners are evenly distributed over the refractory 
walls, 

2. There is no distinction between different elements 
of refoimer tube surface due to geometrical 
disposition or proximity to burners, 

3. The flames are spherical in shape and radiate at 
adiabatic flame temperature, 

4. Burner cup surfaces are no flux surfaces, 

5. All radiations starting from flames reach the 
reformer tubes directly, the only loss to it is 

in passing through absorbing flue g^s media. These 
radiations are completely absorbed on the tube 
surface 

6. The amount of radiations from flue gas and tubes 
absorbed by flames, is negligible. It Aould be 


and 





noted here that the errors incurred due to 
this assmption and those due to assumption 
5 are of opposite signs. 

Based on these assumptions the amount of heat trc 
f erred to reformer tubes by flames can be written as 


1 ?) 

) 


^ t - % -^f ^f ^ 


(2.1 


where lif-g is total number of burners in use, ^ = tc Bg, Dg 
is the diameter of hemispherical burner cup, is the emis 

4 

VI ty of flames, % = is the adiabatic flame 

temperature, is calculated from equation (2,ll) at flame 
temperature, using a mean beam length equal to 0.63 Bg 
i.e. 0.63 times the diameter of spherical gas shape [20], 
Since one element of reformer tube surface is no 
different from another, the transfer of heat to an unit tube 
surface by flames can be represented as 

Ng (l-eg)e^E|^ 


If 


*'t ^ 


(2.17) 


The total transfer of heat to an unit tube surface 
can now be written as 




+ 


If 


t 


N 




t "^t 


( 2 . 18 ) 
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In the above equation it is easy to visualize the 
total transfer as the radiative transfer from flue gas on 
which that from flames have been superimposed, 

(ii) Transfer of Heat from Reformer Tube Inner Wall to the 
Process G-a s: 

It IS assumed that heat transfer from inner tube wall 
is mainly by forced convection. Transfer by radiation and 
conduction have been neglected. It is reported that the 
radial temperature gradient inside a reformer tube is 
insignificant [5] , Therefore, a one dimensional model for 
heat transfer has been used. However, considering the two 
important resistances! one adjacent to the wall and the other 
in catalyst bed, the heat transfer has been calculated from 
the correlation of Seek [ 4 ]. It gives heat transfer coefficien 

h = -f- [2.58 (Re + 0.094(He 

P P P 

(2.19) 

where Re = D G/p, Pr = C D is the diameter of the 

P P P P 

catalyst particles, K, p and C are the conductivity, viscosity 

P 

and specific heat of process gas respectively. Hyman [23] 
has reported that for ring-shaped catalyst particles in a 
reformer tube the heat transfer coefficient is approximately 
40 per cent of the value calculated from equation {2.. 19). In the 
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present irork this is used for all particle shapes and 
the amount of heat transferred from unit surface of inner 
tube wall is 

Ijn = <2.20) 

where is the inner tube wall temperature and T is the 
process gas temperature. 

2.1.3 DESCRIPTION OF THE MATHMATICAL MODEL 

The basis of the mathematical model is a differential 
reformer tube section of length dz filled with nickel catalyst. 
It IS shown in figure 2.2. 

Heat Transfer ; 

It is assumed that there is no transfer of heat in 
the axial direction. This implies that the heat transferred 
to the outer surface of the tube element is conducted through 
it to the inner surface, in the radial direction. The inner 
surface in turn, transfers it to the process gas enclosed 
within the differential section. Since there is no accumulation 
or loss in transfer through the tube element, the heat received 
by it can be equated to the amount conducted through it and 
both of these could be equated to the amount of heat trans- 
ferred to the process gas. These equalities could be 
expressed in terms of heat transfer per unit outer tube 
surface, q, as 
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and 


<1 = 


2 I K (T ~ T ) 
o in 

D in ^ 

0 o' in 


( 2 . 21 ) 


q = D. q /D 
^ in ^in' 0 

where and are the outer and inner tube diameters 
respectively and and are the outside and inside 
surface temperatures. 


Mass and Heat Balance ; 

It is assumed that the performance of a single reformer 
tube IS representative of any other tube in the furnace. The 
temperature, pressure and composition through a general cross 
section of the catalyst bed in a reformer tube has been 
assumed to be uniformly distributed and axial diffusion of 
heat and mass has been neglected. ¥ith these assumptions the 
material and energy balance equations over the differential 
section, dz, (Figure 2.2) can be written as 




d A T P 

1 _ 1 / G 

dz ~ 


G 


for i=l,2 (2.23) 


is defined by 


2 /f 

^ 

g = g.+ N'a.M. 

i=l 


(2. 24) 


and 


^ _ 

dz ~ 41 

1=1 


ah ^ P, 

^ _ t_i.^ _i w 

0 

P 


c ^ 4 in 


G D 


in G 0 _ 




in 


(2.25) 
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where r is the rate of ith reaction, X is catalyst 
bulk density, 0 is the mass velocity, is the heat of 

the ith reaction, G is average specific heat of the process 

ir 

gas, g and g° are the mass fractions of component 3 at any 

0 J 

point and the entrance respectively, M. is the molecular 

D 

weight of component 3 and a is the stoichiometric coefficient 

-L J 

of component j in the ith reaction, with a negative sign when 
the component is a reactant. 


Pressure Drop ; 

Several investigators have proposed correlations for 
finding gas pressure drop for flow in a packed bed. Yen [56] 
evaluated three such methods and found the Brgun [15] 
equation to match his experimental data within 45 per cent. 
The Ergun equation is 


dz 


n ^ n V 


( 2 . 26 ) 


(32.2)(144)(14.7)(3600)^ 

where vis the void fraction and ^ is the density of gas. 

Por hydrocarbon reforming, particle Reynolds number. Re , is 

P 

in thousands. Therefore, the equation being used neglects 
the first term in the bracket. 

It IS noted that the mass velocity, G, is a constant 
throughout the entire tube length and for the moderate pressure 
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and high temperature employed in reformers it can he written 
a s 

a _ M.P. 

V " RT 

where M is the average molecular weight of the process gas. 

The values of and have heen taken from values tabulated 

by G-hoshal et al.[l7]. These have heen expressed as linear 

functions of temperature and used in this work, C has been 

]P 

obtained by averaging specific heats of constituent molecules. 
Heats of reactions, expressed as linear functions of 

temperature have been used [14]. Equations (2.7, 2.8, 2.11, 
2.13, 2.14, 2.18, 2.19 through 2.26) describe the mathematical 
model. 


2.1.4 flAT.GIJLATIQH PROCEDURE 

Process feed hydrocarbon is asstimed to hydrocrack 
completely to methane at the entrance and heat effect asso- 
ciated with the process is neglected. The conversion of 
hydrocarbon (represented as molecule) to methane is 

expressed by the following stoichiometric relation 

3 CH^y + (2-y) H^O (l+y) CH^ + (2-y) CO 

(2.27) 

from which the amounts of CH^ and CO are calculated. Amounts 
of and CO^ in. the ^s are obtained from equilibrium 
considerations of reactions (2.1) and (2.2) respectively. 
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These calculations are repeated till the equilibrium conditions 
of both the reactions are simultaneously satisfied. This 
yields ihe feed gas composition for the primary reformer. 

Adiabatic flame temperature is calculated by balancing 
the sensible heat in incoming fuel and air plus the heat 
generated by combustion with the heat in flames [21], Amount 
of excess air used is calculated from the percent of oxygen 
in the combustion product (flue gas). Measured flue gas 
temperature and calculated adiabatic flame temperature are 
used in equation (2.11) to evaluate emissivities of flue gas 
and flames respectively. In these (Calculations beam length 
for flames is taken to be 0.63 times the diameter of hemi- 
spherical burner cup whereas that for flue gas is 1.8 times 
the gap between refomer tube row and the facing wall [20] » 

To evaluate the heat transferred to a differential 
tube section or the change in the process conditions inside 
it, the tube trail temperatures are to be known. Equations 
(2.18) through (2.22) are used for this purpose. In the 

calculation scheme, a suitable value for T is assumed and 

0 

q IS calculated from equation (2.18). Subsequently, q^^ 
is evaluated from equation (2.22). Equations (2.19) and (2.20) 
in turn lead to the value of Finally, to check the 

validity of assumed value, T^ is calculated from equation 
(2.2l) and in case it differs from the assumed value, the 
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calculated one is used as the new value of for next 
iteratioUc The process is continued till assumed and calculated 
values of agree within the desirable accuracy limits 
( < 0,1°G). 

Inlet process conditions and for the differential 
section being known, equations (2.23) through (2.26) are 
simultaneously integrated using Runge-Kutta method of fourth- 
order [13] to obtain the conditions of outlet which, in turn, 

IS the feed to succeeding differential section. Integration 
IS started, at the inlet end of reformer tubes and continued 
upto the outlet yielding composition, temperature and pressure 
at each point inside the reformer tubes. Values of 

‘If t' ° ‘^in 3-lso evaluated simultaneously 

along the whole length. Size of the integration step 
IS varied according to the length of reformer tubes. 

Calcul ation P rogram ; 

On the basis of the mathematical model described in 
the proceeding section, a calculation program in Fortran 17 
has been prepared for IBM 7044 computer. The calculation 
time for checking the performance of a reformer is 3 minutes 
and 5 secs, when the step size is l/lOOOth of the heated tube 
length. 
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2.1.5 B-SSUITS AND DI3CUSSI0H 

Tlie calcxilation results lor three reformers of different 
sizes, using different hydrocarbon feeds and fuel (naphthas or 
natural gas) and working under various process conditions 
are presented in Tables 2,1, 2,2 and 2,3. Figures 2,3, 2.4 
and 2,5 show the temperature profiles of process gas and 
reformer tube inner and outer wall for different cases. 

Figure 2,6 shows a typical calculated concentration profile 
of the comxDonent molecules along the reformer length (Table 
2,3 Case I), 

The calculated values for process variables are in 
very close agreement with the corresponding measured values. 
Since the rate of reaction is very fast, process gas outlet 
tempera-ture is the most important process variable. The 
maximum difference between the measured and calculated values 
of process gas outleb tem-perature is 3 per cent (Table 2,3, 
Case I). 

The amo;un; of heat calculated as the difference between 
sensible heat of the flames and that of flue gas leaving the 
reforming furnace has been tabulated as total heat absorbed. 
The calculated values of total heat absorbed and total heat 
transferred are in very close agreement. Only in the case 
of short tube reformer (Table 2.2) this difference is 8 per 
cent, in other cases the values are wirhin one per cent. 
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TABLE 2.1 


CALCULATED AND EXPERIIIEMTAL DATA FOR STEAM 


EA.PHTHA PRIMARY 

REFORMER 


Process and fuel liydroca,rbon - Naphtha- [final boiling point 


<210°C, C/H ratio 

Total number of reformer tubes 

160 


Heated tube length, meters 

11.36 


Inside tube diameter, mm 

115.3 


J 

Outside tube diameter, mm 

147.6 


Burner cupt diameter, mm 

356.0 


Catalysb shape and size 

Raschig rings 
I6xl6x6 mm 

Bulk density of catalyst , Kg/m ^ 

1358.0 


Calorific value of fuel, K. cal/Kg 

10200 



Case I 

Case II 

Number of burners in use 

526 

521 

Fuel naphtha flow, kg/hr 

6520 

, 5180 

Flue gas temperature, 

1058 

1010 

Oxygen in flue gas, per cent 

2.0 

4.3 

Total heat absorbed, K, cal/hr 

4. 03x10 

3.406x10 

Process Inlet; 



Process naphtha flow, Kg/Iir 

14100 

11900 

Calculated methane equivalent, 
kg, mole /hr 

695 

587 

Steam flow, Kg /hr 

67200 

56000 

Temperature, °C 

443 

440 

Pressure, kg/cm^ 

24.8 

24.3 
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Table 2.1 (contd) 




Case I Case II 



Plant 

Calculat ed 

Plant 

Calculated 

Process Outlet 





Dry gas flow, Nm^/hr 

61300 

63723 

5410 0 

53850 

Steam flow, kg/hr 

45600 

44188 

35800 

36500 

Temperature, °C 

761 

770.1 

775 

771.5 

Pressure, kg/cm^ 

21.9 

21.75 

21.6 

21.65 

Composition of Gas, Mole per cent 
(Dry basis) 




CO 

9.93 

11.21 

12.6 

11.32 

\ 

CM 

O 

O 

17.43 

16.97 

17.2 

17.01 


64.16 

64.44 

63.3 

64.38 

CH^ 

8.13 

7.38 

6.9 

7.29 

Heat Transfer; 





Heat transfer from 
flue gas, K, cal/hr 


2. 989x10"^ 


2.5025xlo' 

Heat transfer from flames, 
K.cal/br 

1.005x10^ 


0. 8605x10 ' 

Total beat transfer, 
K.cal/br 


3. 9 94x10 


3.563xlo'^ 

Outside Tube ¥all Temnerature,°C 




Reformer Tube Lengtli 

(Practional) 




0^2 

777 

802 

740 

770 

0.4 

832 

842 

825 

813 

0,6 

866 

872 

870 

857 

00 

• 

o 

900 

897 

880 

885 
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TABLE 2.2 

OALCT JLATED AND EXPBRIMEFTAL DATA FOR SHOUT - 
TUBE PRIMART REFORMER 


Process and fuel hydrocarbon - 

Total number of reformer tabes - 
Heated tube length, meters - 

Inside tube diameter, mm — 

Outside tube diameter, mm - 

Catalyst shape and size - 

Bulk density of catalyst ,kg/m^ - 
Calorific value of fuelK. cal/Z@- 
lumber of burners in use - 

Fuel naphtha flow, ICg/hr - 

Flue gas temperature, °C 
Oxygen in flue gas, per cent - 
Total heat absorbed, K, cal/hr - 

Process Inlet ; 

Process naphtha flow, Kg/hr 

Calculated methane equivalent, 
kg. mole /hr - 

Steam flow, Kg/hr - 

Temperature, °C 

Pressure, Kg/cm^ 


Naphtha [final boiling point 
<210® C, C/H ratio 6.27] 

176 

3.28 

90.5 

117.2 

6x6 mm pellets 
1112 
10200 
160 

1814.1 

1045 

3.2 

10.083 X 10^ 

2215.3 

114.72 

13960 

371 

13.7 





TaTole 2. 2 (contd) 


Process Outlet; 

Plant 


Calculated 

Dry gas floWy Efm^/lir 

15800.0 


15236.2 

Steam flow. Kg /hr 

9980.0 


10327.0 

Temperature, °C 

796 


785.3 

Pressure, Kg/cm^ 

13.0 


13.1 

Dry gas composition, 
mole per cont 




CO 

11.15 


10.79 

CO 2 

13.23 


14.36 

=2 

74.00 


72.52 

OH^ 

1.52 


2.33 

Heat Transf e r ; 




Heat transfer from flue gas, K, cal/lir 

- 

6.046 x 10^ 

Heat transfer from flames, E, cal/hr 

- 

3.224 x 10^ 

Total hoat transfer. 

K, cal/hr 


9.27 X 10^ 

Outside Tube ¥all Temoerature , °C; 



Reformer Tube Length 

(Practioiaal) 



0.25 

875 


860.4 

0.50 

904 


897.3 

0.75 

940 


938.4 
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TABLE 2.3 


PLANT AND CALCULATED DATA FOR STEAM-NATURAL 


GAS PRIME! REFORMER 


Process and fuel hydrocarbon 

Natural gas 

[Average composition, mole 
per cent 

CH. - 81.5 


C 2 Hg- 7^0 
C 3 Hg- 5.5 
C 4 H 10 - 4.5 
CO 2 - 1.5] 


Total n-umber of reformer tubes 

200 


Heated tube length, meters 

12.2 


Inside tube diameter, ram 

96.0 


Outside tube diameter, mm 

133.8 


Burner ciixo diameter, imn 

356 


Catalyst shape and size 

- Raschig rings 

16 x 16 x 6 mm 


Bulk density of catalyst, Kg/m^ 

1182 


Calorific value of fuel, K, cal/Nm 

^ 10800 

Case I 

Case II 

Number of burners in use 

584 

568 

Fuel natural gas flow, Nm^/hr 

7786 

706 9 

Flue gas temperature, °C 

1030 

1038 

Oxygen in flue gas, per cent 

2.5 

3.8 

Total heat absorbed, K. cal/hr 

3.96x10^ 

3.69x10' 

Process Inlet; 



Process natural gas flow, Nm^/hr 

13920 

12014 

Calculated methane equivalent, 

Zg. mole/hr 

689.83 

5 98.6 
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Steam flow, Kg/lir 


61000 

55000 

Temperature, °C 


456 

i.4 

447.0 

Pressure, Kg/cm ^ 


35. 

0 

32.3 

PPOCGS?! Olltl ^ 

Case I 

Case 

II 


Plant 

Oalculated 

Plant 

Calculated 

Dry gas flow, ITm^/hr 

62400 

61362.7 

54160 

53342,2 

Steam flow,I{g/h.r 

44000 

44417.0 

39600 

39852.3 

Temperature, °C 

773 

796.8 

793 

784.3 

Pressure, Kg/cm^ 

31.3 

31.56 

28.8 

29.12 

Composition of gas, mole per 
(dry basis) 

cent 



CO 

10.17 

9.56 

9.86 

9.90 

o 

o 

12.27 

12.56 

11.79 

12.01 


67.22 

69. 61 

69.03 

69.19 

CH^ 

9.83 

8.27 

9.32 

8.90 

Heat Transfer: 





Heat transfer from 
flue gas,K, cal/hr 


2.926x10^ 


2.721xlo' 

Heat transfer from 
flames, IC.cal/hr 


1. 004x10"^ 


0. 936x1o' 

Total heat transfer, 
K. cal/hr 


3 . 93x10 


3.657xlo' 

Outside Tube Wall Temnerature 




Reformer Tube Length 

(Fractional) 



0.2 

772 

792.8 

754 

775.3 

0.4 

836 

831.0 

830 

724.2 

0.6 

866 

864.5 

849 

848.0 

00 

a 

891.3 

892.3 

871 

875.8 



000 i 




900 h Outer wall 



Reformer tube length (fractional) 



looa 






^U33j0d (stsoq Ajp) uoipDJi 0101 ^ 


Fig 2-6 -Calculated concentration profile of component molecules m a 



The heat transfer model is more accurately validated 
from the data on reformer tube outside wall temperature, 
at various points along the length, presented in Tables 2.1 
through 2,3 figures 2.3 bhrough 2.5. These figures 
show the calculated temperature profiles of tube wall surface 
and that of process gas along with measured values. The 
agreement is generally very good. However, the measured 
temperature near the ihLot end (0.2 fractional length) is 
lower compared to the calculated value except in the ca.se of 
short tube rofonnoi* (Table 2.2, fig. 2,4). This could be 
because of the assumptions of uniform flue gas temperature 
and uniform disfcribr.tion of burners throughout the wall. In 
actual practice the flue gas temperature at the outlet (refoimer 
tube inlet end) is lower and less number of burner are used 
in this section. This effect is insignificant in the short 
tube reformer because of the smaller length which reduces the 
effect of a localized section. 

In absence of any measurements of composition of process 
gas at various points inside the tube, a proper ;]ustification 
of reaction model is lacking. However, the close agreement 
between all the values measured with the corresponding 
calculated values is a reasonable justification of the 
validity of the reaction rate equations. 

The results obtained do justify all the assumptions 
made to obtain the reaction rate and heat transfer models. 



41 


However, it should he noted that the accuracy achieved in 
the calculations is more "because of the contradictory effects 
of various assumptions made rather than the correctness of 
any single assumption. For example, the number of burners 
used in middle section of the reforming furnace is generally 
more than those in the top or bottom sections. Therefore, 
the assumption of uniform distribution of burners over the 
walls lea.ds to higher values of calculated heat transfer rates 
in the top and bottom sections and lower values in the middle 
section compared to the actual ones, whereas the neglect of 
end effects has just the reverse effect. Similarly the neglect 
of heat loss through the refractory walls and the heat effects 
associated wibh the conversion of higher hydrocarbons to 
methane have opposite heat effects. 
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2.2 SECOITDARY REPOMER 

The feed to secondary reformer consists of PR exit 
to which air is mixed. The SR feed condition is therefore 
dependent on the amount of air used. Too little air means 
insufficient nitrogen and unconverted methane while too 
much air leads to excessive temperature and excess nitrogen. 

The amount of air mixed should be such that at ammonia 
synthesis reactor inlet ratio is 3. Since a small 

per cent of methane remains unconverted, the ratio, 

based on 100 per cent conversion is, taken as 3.1. It is 
further assumed that no oxygen remains after combustion stage. 

2,2.1 PSEP COiroiTIO R 

Mixing of PR exit with air at SR inlet leads to 
adiabatic combustion of 00 and CH^ molecules in the gas. 

The rate as well as equilibrium for the burning (reaction with 
oxygen) of each of them is different and therefore for evaluating 
the exact composition of the product, ail the combustion 
reactions should be considered in simultaneous equilibrium [14], 
However, the rate at idiich hydrogen burns is much faster than 
the other two. Therefore, in this work, it has been assumed 
that the extent of combustion of 00 and CH^ is negligible 
compared to that of Close agreement was found between 

the combustion product temperature calculated on this basis 
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and that calculated by Davis and Lihou [14]. this 
model the amount of air required is calculated as follows; 


air 



(2.28) 


where ^CH inlet rates of air, CO, 

H 2 and CH^ respectively, Nm^A^^ and and are the mole 

fractions of N 2 ^2 air. Amount of H 2 , 1^2 and Ar in 


the combustion product are obtained from the following 


relations 



(2.29) 

(2.30) 

(2.31) 


where and are the rates of flow of H 2 , ^2 and 

Ar in the process gas after combustion is complete. Amounts 
of other components in the PE ©zlt remain unaffected since 
they do not participate in the reaction with air. 

Temperature of the process gas after combustion, and 
before it reaches SR catalyst bed, is calculated by equating 
the sum of sensible heats in air and PR exit and the heat from 
combustion to the sensible heat in combustion product. This 
IS calculated by a trial and error procedure [21]. 
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2-2.2 RATE EQIJATIQITS 

As in PR, a nickel catalyst is used in this reactor 
also with the only difference that this catalyst is capable 
of standing relatively higher temperatures. The pair of 
equations used to describe the reactions is same as that 
for PR catalyst i.e. 


CH^ + H^O — 


CO + H^O 


CO + 3H, 


CO2+ 


( 2 . 1 ) 

(2.2) 


The same kinetic expressions (Equations (2.7) and 
(2.8 )) can also be used here. However, at the SR temperature 
condition, fate of reaction (2.1) is so fast that it could 
be assumed that the methane molecules hydrocrack instantaneously 
at the catalyst surface. The surface concentration of methane 
IS, therefore, taken as that in equilibrium corresponding to 
catalyst surface temperature. 

In such situations the rate of reaction is governed 
by the rate at which the reactant diffuses to catalyst surface 
from the bulk gas [54] which could be expressed as follows; 

where C concentrations in the bulk gas 

CH^ ^^4 

and at the surface respectively, Eg.mole/Nm^; is the mass 
transfer coefficient between bulk gas and solid surface, 


meters/hr. and is the external surface area per unit mass 
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1 

of the catalyst pellets, m^/Zg, of catalyst. C^-o- and the 

4 

catalj'-st surface temperature, are evaluated from a.n energy 
■balance oxi the pellet which could "be written, after proper 
manipulations and approximations [34] as follows 

T-T® = 0.7 X (- X -i (2,33) 

C . 

where C is the total concentration of the process gas, Eg,mole/m^. 

Tor the solution of a'bove correlation, which involves trial 

and error solution, reaction (2.2) is assumed to "be absent. 

Rate of shift reaction is described by equation (2,8) 

i . e . Y • T 

“^00 

r = exp(- + 8.23)fP (Iqq" 

^ 0 ( 3.2 -^ 2 ^ 

2.2*3 DESCRIP TION OR hATHEIlATI CAD I-IQDEL 

The basis of the mathematical model i.s a differential 

section of length dz, of the adiabatic catalyst bed. The 

assumptions and material, energy and pressm’e balance equations 

over this section are same as 'chose for PR catalyst sec cions 

except that heat transfer term in equation (2.25) is zero. 

Equation (2.25) can be rewritten for this case as 

y ("^) (2.34) 

dz ^ Q ' G 

1=1 p 

Equations (2,8) , (2. 23) , (2.24), (2,26), (2,32), (2,33) and (2.34) 
adequately describe the perfo'cmance of SR adiabatic catalyst 


bed. 
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Simultaneous solution of equations (2.25), (2. 24) ,( 2. 26) 
and (2.34), using Runge-Zutta method of fourth order [13], 
yields temperature, composition and pressure bhroughout the 
catalyst bed. 


2.2,4 RESU LT S ARR DISCUSSION 

The calculation results for three reformers are 
presented in Tables 2.4, 2.5 and 2.6. figure 2.7 shows the 
variation of calculated temperature along the catalyst bed 
length alongwith the measured values at the inlet and exit 
(Tabic 2.4, Case l). 

The calculated SR inlet temperature (after mixing vith 
air) agrees well with those reported by Davies and Lihou [14] 
who have considered simultaneous combustion of H„, CO and CH,* 

It is to bo noted that in the present calculations, combustion 
of only is considered being the fastest among the uhree. 

The calculated values also agree with the measured SR inlet 
temperatures which are available (Table 2.4). These agreements 
validate the assumptions made. 

The agroemont between measured and calculated composition 
and temperaturo for all cases, at the SR outlet, shows the 
correctness of tho assumptions made in tho calculation of 
steam-methane reaction over SR catalyst. However, the validity 
of the method of calculation cannon be ascertained because of 
lack of data on temperature or compositioi at intermediate 
points in the SR catalyst bed (Figure 2,7). 
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Fig. 2-7 -Temperature variation in 5R catalyst bed. [Table 2-4,CaseI] 
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means low carbon monoxide consent in the exit gas. 

Some of the old units have only high temperature shift 
reactors. Multiple beds with interstage cooling are used -with 
progressively lower temperatures in the subsequent beds. 
However, the lowest temperature has to be kept high enough 
in order that the HT catalyst is still effective. Current 
practice is to provide an LT reactor in conjunction with a 
single bed HT reactor. Use of IT reactors has several 
advantages [16]. Yery low (<[0.2 per cent) concentration of 
carbon monoxide in the exit from LT can be attained. With 
such low concentration of CO, methanation catalyst can conve- 
niently be used for removal of last traces of carbon monoxide 
and dioxide without involving high purge loss. 
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IHCtH TMPE PlATURE VATER-gAS SHIFT REACTOR (HT) 

For resc'bion over HT catalyst (chromia promoted iron 
oxide) some workers propose a first order [3,26,28] or a 
second order [29-31] rate equation to explain the conversion 
tut many others [4-6,8, 9,15, 20, 22-25, 32, 37] propose 
different rate equations. Each of the latter rate equations 
IS of comp] ex order and different from those proposed hy 
others. It l-'o experimentally been observed that catalysts 
of different manufacturers give differenx rate equations even 
under identical conditions [15]. The specific surface area 
of the catalyst is reduced considerably during use due to 
sintering [18] and the surface composition of the catalyst 
changes depending on the composition of the reacting gases 
which are a mixture of oxidizing and reducing gases. All 
these may help in qualitatively explaining the difference in 
rate equations obtained by different workers but it is not 
veiy viseful m getting a general rate equation which is not 
available in literature. So in this work the approach is to 
start with a rate equation which explains the conversion under 
some constraints. Each of these constraints are subsequently 
eliminated, one at a time, by introducing suitable correction 
terms. These are mostly empirical correlations obtained from 
experimental results of different workers on the same catalyst. 
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3.1.1 DBYELOPI IMT OB T HE BATE EniJATIOI 

The rate equation most frequently referred to in 
literature is first-order or pseudo-first order in carbon 
monoxide concentration or partial pressure [3,26,28] an(^ all 
attempts to correlate the rate data by power function models 
give reaction rates which are almost consistently first order 
with respect to partial pressures of carbon monoxide [5,8,15]. 

So a first order rate equation of the form given by equation 
(3.2) has been selected. 

r = h’ (Eqq- Eqq) (3.2) 

where, r = rate of reaction, cc/hr.gm. of catalyst, 

k' = rate constant, cc/hr. gm. of catalyst, atms, , 

= partial pressure of CO, atms. and 
P*Q = partial pressure of CO in equilibrium conditions, ati 
It IS noted that this rate equation has been used under low 
pressures (of the order of 1 atmosphere), high sfeam to carbon 
monoxide ratio and below a temperature of 500°C. In commercial 
reactors latter two conditions are met but the operating 
pressure is generally high (l5 to 30 atms). So the above rate 
equation is used at a constant pressure of 1 atmosphere, and 
the variation in pressure is eccotintod for separately. At one 
atmosphere pressure, equation (3.2) can be written as 


r = k (xgo - x*q) 


(3.3) 
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where, k = rate constant, cc/hr, gm. of catalyst. 


Xqq = mole fraction of CO, 


= mole fraction of CO in equilibrium conditions 


xJt X X* 


"H. 


CO. 


* T7- 

x„ ^ X K 

H_0 eq 


and = equilibrium constant 

The rate constrnt, k is generally expressed in the 
form of Arrhenius Equation as 


k = a e 


-E/RgT 


(5.4) 


where, a = pre-exponential factor, cc/hr. gm. of catalyst 

E = activation energy, cals/gm.mole 
R = universal gas constant, cals/gm. mole and 

o 

T = absolute temperature, °K 
The value of activation energy, E, to be used in this 
equation is reported from 13,000 to 32,000 cal/gm mole [12, 15, 
21,30] but mostly it lies in the range 21.4 to 27.3 Kcal/gm 
mole [35]. Eor the present catalyst (of P and D, FCI) the 
value of E is calculated as follows: 

A pilot plant study using pellets of this catalyst 
and feed gas containing high CO content (40,6 per cent, dry 
basis) IS reported in literature [13]. In this, the calculated 
value of energy of activation for the catalyst without taking 
into account any mass or heat transfer resistance between the 



"bulk gas and catalyst' pores i.e. the apparent energy of 
activation for the catalyst pellets, is 

E’ = 15.88 Zcal/gm mole (3i5) 

For commercial shift reactor operating conditions, 
heat transfer and the external mass transfer resistances may 
be considered insignificant but the intra -pellet mass transfer 
resistance is very high [28]. In such a case, the true energy 
of activation, E, can be taken as twice the apparent energy 
of activation, E' [36J. Therefore, 

E = 2B’ = 27.76 Zcal/gm mole (3.6) 

In the above evaluation any temperature difference 
between the catalyst pellets and the bulk gas or between 
different points in a catalyst pellet has been taken to be 
insignificant. The close agreement in value of E thus obtained 
and those reported by different workers [14,21,35] who used 
small [ <1^ 1 mm size] catalyst particles, lower pressures and 
lower CO concentrations in their studies, suggests that it is 
so. However, if some temperature difference exists, its 
effect IS absorbed in the rate equation and hence need not be 
considered separately. 

The pre-exponential factor, a, is not a constant. It 
depends on temperature because the specific surface area 
varies with temperature [18], In thiB work, a has been 



evaluated at a temperature of 350°C by equating experimentally 
determined intrinsic rate constant [33] to that given by 
equations (3.4) and (3.6). This value of 

a = 2.865 X 10^^ cc/hr gm. of catalyst (3.7) 
has been used throughout this work. * 

Error incurred in the calculation of a catalyst pellet 
activity due to taking a constant value of a is corrected for 
temperature as discussed in the section dealing with catalyst 
age, Erom equations (3.3), (3.4), (3.6) and (3.7) 
the intrinsic rate of reaction for the catalyst is given by, 

r = 2.865 X 10^^ exp 2 : (x^^-xj^) 

(3.8) 

Rate of reaction for the pellet is obtained by incor- 
porating the effects of diffusional resistances, age of the 
catalyst, pressure and H^S concentration in the reacting gases. 
These factors have been discussed in the above order in the 
following paragraphs. 

( a ) Diffusional Resistances ; 

As mentioned earlier only intra -pellet diffusional 
resistance is to be considered. This can be done by using a 
suitable effectiveness factor. Eor calculating effectiveness 
factor, catalyst pellets have been considered to be isothermal[lC 
and Wheeler's method [42] has been used. This gives 




Effectiveness factor. 


•where, 


Eff = 

(3.9) 

h = , 

(3.10) 

0 e 


T 


L = f2 

^x 

(3.11) 

“e- +V7 Bb 

(3.12) 


7p= volume of the catalyst /pellet , 
S = external surface area /pellet, 
r^ = pore radius, 

D = effective diffusivity, 

Bg = bulk diffusivity and 
= Knud sen diffusivity 


It has been sho 'wn [28,33] that for the shift catalysts 
at low pressure, Knudsen diffusion is the controlling mode of 
diffusion. So a method based on pore size distribution has 
been used. This method as well as the pore size distrihution 
for the catalyst is reported elsewhere [33]. In this method 
effectiveness factor, Eff(r^) corresponding to different 
pore sizes (^q) is calculated. This is multiplied with the 
total surface for the pore size, A^. This gLves the effective 


surface area 



L 


corresponding to pore size, r^. is calculated for values 

of r^ in the range of 90 to 250® A, •wh.ich accounts for ninety 
per cent of the total surface area. Now effectiveness factor 
IS calculated at one atmosphere from the following relation 


Eff = 


To tal effective surface area 


Total pore area 
250°A 


r =go 
o 


(3,13) 


250°A 

A 

^ P 
r =90 
0 

(h) Age of the Catalyst ; 

A catalyst on-stream continuously lo£«-s. its activity 
with time at a rate winch depends on the operating conditions 
to which it is subjected. Above 500®C this loss in activity 
in the case of a shift catalyst corresponds to the loss in 
specific surface area, A. A is expressed in terms of the 
operating temperature, T, and the age of the catalyst, t, as [18] 


A' 


n 


.n 


= Kg T e 


(3.14) 


where, A^ = specific surface area of a fresh catalyst. 

Kg = a constant and 

Eg = activation energy of sintering. 

Eor large values of r, this is reduced to a simpler 


form: 



Do 


A = X X exp^®s/\^^ (3.15) 

The above equations show that the surface area for a 
given process time, t, decreases with increase in temperature 
and for large process time (t) (i.e. when eqn.5»15 holds) 
the variation in area is small for large variations in t; 
since n lies between 6.2 to 7.3. But the shift -convert ers are 
operated below 500° 0 where the correspondence of A to catalyst 
activity does not hold. The actual loss in catalyst activity 
IS much more than the loss of surface area [13,28], This 
has also been pointed out in the original work [18] on reduction 
in surface area due to sintering. Therefore, an empirical 
correlation has been obtained, which is based on a study of 
the quantitative variation in the activity of pellets of this 
catalyst with temperature and time [13], A typical curve 
of catalyst activity vs time, from this study, at any parti- 
cular temperature shows (Figure 3,1) that the activity falls 
very sharply in the initial stages. Then almost abruptly, 
the rate of fall in activity is reduced to a very low value. 

The activity of the catalyst at the turning point M has been 
compared with those calculated from intrinsic rate constant [33] 
for the catalyst reduced at the temperature for which the 
curve PMl containing M has been obtained. Both these values 
turn out to be equal within the range of experimental accuracy. 
This implies that point M's on activity curves for different 
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Fig 3 1 - A typical curve of variation in catalyst activity 
with time at constant temperature 



temperatures represent the activities of the catalyst 
completely reduced at respective temperatures, Erro,c3.ta?.yst acti- 
vity is considered onlj?- after it has completely "been reduced 
1, e, the sharp fall period is not considered. There is 
not much loss of generality due to this since the catalyst at 
test takes a few days to stea.dy its rate (i. e. to reach M) 
whereas the total estimated life of the catalyst is three 
years. 

As the activity of a reduced catal3rst is a function of 
the reduction temperature, a co'^relation of the relative 
activities of the catalyst nellets reduced at different 
temperatures is required. Relative activity (Rg^(T)) has been 
defined as - - 


EJT) 


Activity of a catalyst pellet reduced and working 

at T^K 

Activity of the same catalyst reduced at a standard 
temperature T and working at T°K 

^ ( 3 . 16 ) 


Taking standard temperature T^=:623.2°K (350°0), 


£ (rp) _ Activity represented by M for T^K 

^ ' Iff z 2.865 X iC^^x exp( ~-^j — ) (3.17) 

e 

This relation gives R„(T) for different values of T. 

3 , , 

To correlate R (T) with T, an exponential decrease in catalyst 

St 

activity with increase in reduction temperature (as in Bqn.3.15) 
has been assumed as a first approximation i.e. 



[0(1/T-1/Ig)] 

(3.18) 

A plot of -In vs T gives a straight line. This 
validates the assiomption and gives 

R^(T) = exp (-10.831 + 6750/T) (3.19) 

o. 

It is to be noted that one needs to consider the 
activities of only those catalyst pellets which are reduced 
at Operating temperature^sxnce one cannot do “otherwise but to 
assTome that in a converter the catalyst pellets are reduced 
at the operating temperatures only. In the present work the 
possibility of it being subjected to a higher temperature» for 
a period long enough to reduce it (even partially) at that 
temperature is not considered. 

Now one can find out the activity of a reduced catalyst 
pellet at atmospheric pressure and any temperature (i.e. M) 
from Equations (3.8), (3*13) and (3.19). Any further fall in 
catalyst activity (i.e. along ML) with respect to T and t can 
be expressed in terms of another factor (aging factor) which 
has been correlated as 

log^Q Ag^ = (14.66 X 10”'^ -2 X 10“^ x T) x t (3.20) 
from data obtained from the same study [13] referred above. 
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( c ) Pressure ; 

All investigators [5,8,12,29] "but one [38] have 
reported that the rate of reaction increases with increase 
in pressure. In an investigation over the catalyst [12] 
the ratio by which the rate increases has been related to 
pressure as 

Pressure factor, P^^ = “ P/250) (3.2l) 

This IS in voxy close agreement with results of 
others [3,29]. 

(d) H^S Concentration: 

The activity of an iron-oxide catalyst (HT . shift 
catalyst) is reduced by the presence of H^S in the reacting 
gases due to reversible poisioning of the catalyst. The extent 
of the reduction in catalyst activity depends on the concen- 
tration of H^S in the reacting gases. Many workers [7,8,28,38] 
have studied this effect quantitatively. The one reported for 
this catalyst [7] has been used to correlate the reduction in 
activity with H 2 S concentration (in ppm) as 

fg = -0.276 log3_Q([H2S] + 2.78) + 1.127 

(3.22) 

where fg stands for the factor by which the rate is reduced 
and [H 2 S] represents concentration of [H 2 S] in ppm. 



Now the final rate equation can be whit ten as, 




r = Eff X 2.865 x 10^^ x exp (~^^) ^ Eg_ x x P^' ‘X 
* § 

^ ^^C0“^C0^ (3.25) 

which gives the rate of reaction over the pellets of the 
catalyst, whose specifications are given below, under all 
operating conditions which_may exist in "aBy' commercial H3? 
. shift cohve'rter. 


Catalyst snecification Fresh 

Chemical composition 

^^ 2^3 cent 

Cr^Oj 8 per cent 

Real density, gm/cc 
Pcnsity of the pellet, gm/cc 
Specific surface area, m^/gm 47.4 


Reduced at 350^C 


4.561 

2.225 

31.2 


Porosity ; 

^200 °A 0.2436 cm^/gm 

200 '^A-75,000°A 0.0818 cm^/gm 

;>7500 °A 0.0019 cmVgia 

Total porosity 0.2901 cm^/gm 0.5273cm^/gm 


3.1.2 DESCRIPTION OF THE MATHIMATICAl M ODEL 

The following mathematical model ha,s been developed 
for a high-temperature shift reactor which may have single or 
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multiple adiabatic beds of catalyst with, interstage cooling. 

One such reactor having three adiabatic beds with interstage 
cooling IS shown in Figure 3.2. 

Unifom distributionof temperature and composition 
throughout a general cross-section of the catalyst bed has 
been assuned and axial diffusion of mass and hoat has been 
neglected. ¥ith these assumptions the material and energy 
balance equations which describe th«^ composition and temperature 
of the reaction system along an adiabatic catalytic layer can 
be written as 


r» 

dv ~ & 


(3.24) 


dv 





G 


for j=l, 


4 


(3.25) 

(3.26) 


where 


V 

G 

r ’ 


X . ,x^ 

3 3 


oc. 


= volume of catalyst 
= volume flow rate 
= reaction rate, 

= heat of reaction, 

= average specific heat of reacting system 
= mole fraction of component j at any point and 
the entrance respectively, and 
= stoichiometric coefficient of 3 , with a negative 
sign when it is a reactant. 



(unconverted) 


■1st catalytic bed 


2nd catalytic bed 


Sh 1' 

I l ji 


3rd 

catalytic bed 


•Inter stage 
cooling heat 
exchanger 


Converted qas 


2 -'Triple adiabatic bed reactor with interstage coolin 



Since the reaction takes place at a moderately high 
pressure, reacting gaseous system can not be considered ideal. 
So the pressure dependence term has been included in the 
expression for specific heat [19] of each component using 
data [17] on different pressures (l-40 atms). The standard 
heat of reaction at 25*^0 and 1 atm pressure has been 
calculated from the heat of formation data [19] which has 
jjhen been generalized for calculations at hi^er temperatures 
and pressures [17]. 

The composition and temperature of the reacting system 
at the inlet of any layer being known, the material balance 
equation (3.24) is solved to obtain the extent of reaction 
which in turn gives the mole fraction of individual components 
(eqn. 3.26). Solution of the energy balance equation (3.25) 
gives the outlet temperature. Runge Kutta method of fourth 
order [11] has been used for the simultaneous solution of mass 
and energy balance equations. The size of the integration 
step has been varied according to the volume of the catalyst 
used. 

Calculation Program ; 

On the basis of the mathematical model described in 
the previous paragraph, a calculation program in Fortran IV 
has been prepared for IBM 7044 computer. The calculation time 
for checking the performance of a triple~bed reactor with 



1 n't 01 * 81:5 go cooling 3.s 53 seconds vlieii tlio step size for eacli 
bed is l/lOOtli of tte bed voluric. Tins xorogran lias been used 
only Tor clicclciiig tl e perTomance of sinft reactors but it 
can conver lently be used for desijr purposes as well, 

3.1.3 ?u S, SUITS AFl D I- SCUSSI OI'T 

The calcul8vtio:i results for five different units having 
single or jaulfciple beds and t^ori-ing iindei different conditions 
of tonpcrabu-c (603 to pressurr^ (l6.4 to 30 Kg/cm^), 

compos3- tion of' tie vea.cbing s 3 ^s+?pi (CO ■“ 9.0 bo 40,6 per cent), 
ca'Gal 3 -’'ot -i'- e (oO bo 560 days), H 2 S concencrstion (less than 
0,5 bo 50 p-'ui) axid steam bo 00 ratio (5 to 21.7) are presented 
tilong m’ bhe plant data in Tables 3.3. to 3.4 and figures 3.3 
CO 3.5. 

For t^e purpose of discussion vhe reactors are 
considered separe cely (based on number of catalyst beds), 
Trinle Bed Rea ctor; Tlin "’arimum difference in concentration 
of 00 calcul"'ted 3 ron that rieacured exicts am the outlet from 
1st bed (Tabic 3.1, Fig. 3.o). Based on the model, the volume 
of catalyse in Ist bed requiz’od for the measured conversion 
( 14.2 per cent of Co, drj*- basis) is 7 per cent less than the 
actual voIujUG. denerally the experimental error upto + 5 per 
cent is allowed. The resiat is slightly beyond that limit. 

But this would be of much loss st yeuiicance if the error in 
double bed reactor, to b-^ iircussed next, is considered. 
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Tlie maximum difference in temperatnre is 11,1°C (Talkie 3.1, 

Fig, 3. 3) which IS only 2.4per cent. 

Double Bed Reactor ; Calculated volume of catalyst in the first 

bed required to obtain the desired conversion is 4 per cent 

more than the actual volume. This error compared to that in 

the Triple bed reactor makes it clear that the error is 

rathei* random. Temperature calculated and measured differ 
o , . 

only by 1.7 C (Table 3.2) which is well within the experimental 
error. 

Singl e Bed Reactors ; Although the difference in calculated 
and measured composition and temp ere tiire at the exit is not 
high (Table 3.3), this alone is not sufficient, particularly 
when the equilibrium is closely approached at the exit. However, 
it can easily be shown, from thermodynamic considerations, that 
if the temperature at any section of the reactor is kno>na, it 
can be used ro predict the composibion of the reacting gases. 
Hence, following the temperature profile is equivalent to 
knowing the conversion profile. This is particularly important 
in the absence of any method of knowing the concentration 
profile. In Table 3.4 the calciilated and measured values of 
temperatures at different sections of the single bed reactors 
are presented. Pignre 3.5 shoT7S the temperature profile of 
Reactor I (in Table 3.3). The curve shows that the measured 
and calculated values are in very good agreement, within 2 per 
cent, and randomly distributed. 
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The results very clearly show that in general, the 
calculated values are in very good agreaaent with the plant 
data. This confirms that the developed rate equarion as 
well as the assumptions made are valid under all plant 
conditions and effect of diffusion of reactants and products 
inside the catalyst pores has been correctly accounted for. 



3.2 L0¥ TEMPBEIATURB WATER-GAS SHIFT REACTOR 


The studies pertaining to the preparative and kinetic 
aspects of low temperature shift -catalyst ( essentially a 
copper-zinc oxide system with various promoting components) 
operating in the range 180 -250°C are relatively recent [14, 16, 
21, 39,40,41] and not as exhaustive as on HT catalyst. The 
kinetics of the reaction is reported to be same as that over 
a HT catalyst [14]. Kinetic studies leading to rate equations 
different from that for HT catalyst are also reported [39,41]. 
However, in general, a 3?ate equation similar to that for 
reaction over high temperature catalyst is assumed to evaluate 
the constants and t o compare activities of the two catalysts 
[1,14,21]. These are used in all calculations for LT 
reactors [27]. 

Equation (3.3) used to describe the rate of reaction 
at atmospheric pressure over a HT catalyst is used in this 
case also. To calculate the effective rate over catalyst 
pellets, the effects of intra particle diffusional resistance, 
age of catalyst and pressure are separately accounted for 
as in the case of HT catalyst. 

Thus, the intrinsic rate at atmospheric pressure for the 
low temperature catalyst is expressed as 

r = - x*o) 

k = a 


where 


(3.3) 

(3.4) 






The value of activation energy , E, is reported between 
20^600 to 26,800 k. cal/kg.mole [1,14,21] for the catalyst 
over •v^ich the model has been tested, the reported [1] values 
are 

E = 20,960 K. cal/kg.mole 

and a = 2,955 x 10^^ cc of CO/hr,^ of catalyst 

at atmospheric pressure. 


With these values, the intrinsic rate at atmospheric 

pressure for the catalyst can be represented as 

2 0.960 

r = 2. 955 X 10^^ e' (x^^q-xJq) (3.27) 

The specifications of the catalyst ov-er which the rate 
IS described by equation (3.27) are as follows: 


Catalyst Snecifl cations ; 

Chemical composition 
CuO 
ZnO 

Real density, gm/cc 
Density of the pellet. 
Specific surface area, 
Porosi tv 
<200 

Total porosity 


- 33 per cent 

- 66 per cent 
5.09 

gm/cc 2.575 

m^/gm 27*6 

0. 208cm^/gm 
0. 243cm^/gm 
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( a) Piffusional Resistance ; 

Knudsen diffusion is the controlling mode of intra- 
particle diffusion. Therefore the effectiveness factor is 
calculated by the method described earlier using equations 
(3.9) through (3.13). 


( b ) Age of the Catalyst ; 

Since the reaction over this catalyst takes place at 
a low temperature, there is no loss due to sintering as in 
the case of HT catalyst. The ;j.oss in activity due to age is 
accounted for by the aging factor 



^q(4.66x10”‘^-1.6x10“^T)x t 


(5.28) 


( c ) Pressure ; 

Effect of pressure on the rate of reaction has been 
studied by some workers [1,39]. The one based on the 
catalyst under study is however inconclusive andhence on 
the basis of similarity the pressure effect is accounted 
for by the factor obtained for high temperature catalyst 
1. e. 

Pj = P^°-^ ~ T/250) (3.21) 

The agreement of plant data and calculated results 
justifies its use. 

Now the final rate equation can be expressed as 



, ^ - 20960 

r = Eff x2,S55 X e R T ^ 

s 


^gf ^ ^ (^GO ^00^ 


(3.29) 


De s cription of Mathematical Model; 

Assimaptions made in the case of HT reactor are valid 
in this case also and the same set of equations (3*24 - 3.26) 
are used to describe the reaction system along an adiabatic 
catalytic layer. Heat of reaction is calculated from the 
same correlation and integration is also performed according 
to the same scliome. 

HBSHLTS add DISGU SSIOF 

The calculation results for three different reactors 
using catalysts of different age (60 to 560 days) and -working 
under different conditions of pressure (l6.0 to 29.3 atms) 
and temperature (183 to 220°C).are presented along with the 
plant data in Tables 3.5 and 3.6, Pigure 3.6 shows that 
the calcvilated temperature profile in Reactor II (Tables 3.5 
and 3.6) along -with the plant data. 

The results show a very good agreement between the plant 
data and calculated values, the maximum difference being less 
than 2 per cent. However, the composition and temperature 
changes involved in these reactors are so small that the 
correctness of the model is not realized very easily as in the 


case of HT reactor. 
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lOMENGLATURE 


4 


A 


\ 

E 

E* 

Eg 

G 

AE 

d 

k 

k' 

P 

P 

R 

S 

r 

r ’ 

T 


Specific surface of catalyst, m^/gm 
Specific surface of fresh catalyst, m^/gm 
Frequency factor, cc/hr, gm. of catalyst 
Specific heat capacity of reacting gaseous mixture, 

Cals/gm mole °K. 

Effective diffusivity, cm^/sec 
Bulk diffusivity, cm^/sec 
Knudsen diffusivity, cm ^ /sec 

Intrinsic energy of activation, cals/gm mole °K 
Apparent intrinsic energy of activation, cals/gm mole °K 
Energy of activation for sintering, cals/gm of catalyst, °E; 
Yolume flow rate, I'm ^ /hr 

Heat of reaction, cals/gm mole of CO reacted 

H'um'ber used to represent different components 

Hate constant, cc/hr. gm. of catalyst 

Rate constant cc/hr. gm of cataljrst, atms 

Total pressure, atms 

Partial pressure, atms 

Gas constant, cal/gm mole 

Rate of reaction, cc/hr gm of catalyst 

Rate of reaction, cc/hr cc of catalyst 

Pore radius, cms 

Temperature, 
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V Voliune of catalyst, cc 

X Mole fraction 

C,K„,n Constants 
s 

a Stoichiometric coefficient (-ve value is used for 

reactants and +ve for products of the reaction) 

X Age of the catalyst, days 

SuhscriTjts 
0 Component o 

CO Component carbon monoxide 

Sunerscrint 
0 

* 


Inlet conditions 
Equilibrium conditions 
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CHAPTER 4 


SIMULATION OE AMHOEIA SYNTHESIS REACTORS 

Tlie exit gas from methanator contains nitrogen and 
hydrogen in appropriate proportion for the ammonia synthesis 
reaction. This gas also contains some methane and argon 
as aiKayb components. The gas mixture is compressed to 
200—550 atmospheres and fed to ammonia reactor along vath the 
recycle gas. The reaction "between nitrogen and hydrogen is 
one of the simplest kinetic reactions. The synthesis is 
straight forward, there is no side reaction and the product 
is stable* The physical and thermodynamic properties of the 
reactants and products are well known [15]. However, the 
mechanism of this reaction over the synthesis catalyst [Iron 
Catalyst] is not well understood. This has led to numerous rate 
equations, all of which are of complex order. Because of 
complexity of these rate equations it is difficult to account 
for the diffusional resistances to the transport of reactants 
and product in catalyst pores. This, in addition to limited 
reliability of rate equation, makes it difficult to have a 
mathematical description of the processes taking place inside 
an ammonia synthesis reactor, 

Maoor changes have taken place in the design of ammonia 
synthesis reactors since the first commercial production 




started in 1925 [19]. Most of these changes have been based 
on historical plant data rather than an insight into the 
physical and chemical processes taking place in the reactors. 
However, the nse of computers in design, optimization and 
control made it necessary to have a mathematical description 
of the process. Simulation models for ammonia synthesis 
converters of different types have been developed for design, 
optimization [26, 35] and control [31,32] purposes. 

To describe the reactor operating conditions as accu- 
rately as possible, the simulation model snould take into 
consideration all t]ae physical and chemical processes taking 

* ' I 

I 

place in the reactor. In order to avoid the complexity 
resulting from such a consideration, earlier workers have 
attempted only approximate simulations. 

The mathematical model described here considers all 
physical and chemical processes in the reactor. Method to 
solve transport equations to evaluate effectiveness factor, 
upto desired accuracy is developed ‘and used in the model 
calculations. 

rate exp r essions 

The literature contains innumerable rate expressions. 

Those reported till early thirties have been summarized by 
Frankenburg [14] and Emmett [11, 12]. In 1940 Temkin and 
Pyzhev [37] developed a rate equation which offered a satisfactory 
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kinetic approach to the synthesis and decomposition of ammonia 
over doubly promoted iron catalysts. Since then this rate 
equation as such or in modified foims has been most extensively 
used, although some doubts about the generality of the equation 


have been raised [1, 13,17]. The modified form of the Temkin 
equation [10], used in this work, is as follows: 




_M^^l-aj 
-p 3 


(4.1) 


where r^jg = reaction rate, kg. mole, ffi^/hr/m^ of catalyst, 
K 2 =veljocity constant of the reverse reaction, kg. mole/hr/m^, 
f , f jj » fjjjj = fugacities of nitrogen, hydrogen and ammonia 
respectively, k^^ = equilibrium constant of the reaction: 1.5 H 2 
t 0,5 and a=constant. According to some workers 

a=0.5 for all iron catalysts [25,36,37] whereas others obtained 
values ranging from 0.4 to 0,8 [2-7,10, 13, 16, 18, 22, 23, 24, 
28, 29, 33, 34]. 

Having found different values of a, many authors suggest 
that it depends upon catalyst characteristics. It appears 
reasonable, considering that catalysts of different make differ 


in their promoter contents and physical characteristics. 

However, difference of the same order in the value of a has 
been reported for the same catalyst operating under two different 
conditions (before and after thermo re si st an cy test). Even for 
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two arbitrary sets of runs over the same catalyst, values of 
a show similar differences [16], 

Oenerally to evaluate a, is expressed as a function 
of temperature in Arrhenius equation form i,e. 


^2 = 


K 


20 


-Ej/E^T 


(4.2) 


where is a constant, is the energy of activation, 

R IS gas constant and T is the absolute temperature. The 
S ^ 

experimental results are then fitted with E2 a-nd a as 

parameters. Such a method of evaluation could lead to 
variations in values of a with a corresponding change in 
values of E2 and E2Q. It has been suggested that there is 
a linear relationship between the value of log ^2 
This shows that to differentiate the effect of catalyst chara- 
cteristics or promoter content from those associated with 
measurement .errors or method of calculation is difficult. 

Based on this it has been suggested that a constant value of a, 
independent of catalyst make, could be used [16]. However, 
for a fixed value of a, other parameters K20 
available for the two catalysts over which the present model 
has been tested. The values used for Montecatini Edison 
catalyst are a = 0.647, B 2 = 38690 Real /kg mold, log K 2 o= 14 » 206 
and the corresponding rate equation is 
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*1 


^NH ~ (2.305 X 14 . 206 - x 

3 § 

f ^ 

- ( —3) 
f/t 

2 


IK ^ ^ (-T^) 


■N 


0.647 




M^^0.353 

3 

IL 


(4.3) 


Eor Haider Topsoe catalyst, a = 0.692, $2 = 42893 and 

>• i I 

log K2Q = 15.2059 and the rate equation is 

(2.303 X 15.2059 - ^ 

f ^ 


X (-|a )0.692 _ (l!i)0.308^ 


'EH. 


■?2 


(4.4) 


All the above values of a, ^2 and log K 2 Q have been 
obtained from the published work of Guacci et al. [16], 


4.2 EFPECTIVMBgS PACTOR 

The above expressions represent the intrinsic rate of 
reaction i.e., the rate of reaction for small particles in 
which there is no resistance to transfer of mass or heat to 
the active surface. The large size industrial catalyst 
particles (6-12 mm for axial flow converters) are, however, 
subject to diffusion restriction in their pore structure. 

This effect can be taken care of by effectiveness 
factor, r]f which is defined as the rate at which the reaction 
occurs in a pellet divided by the rate at which the reaction 
would occur if the concentration and temperature throughout 
the pellet were same as those at the outer surface i.e. 



n 


[Surface area of 
pellet] 


[Molar fliix of 
component i 

surface! ^ 

[Volume of pellet] [Rate of formation of 

component i at surface comp. , 

(4.5) 


Several authors [4,10,20,21,28] have calculated this factor 
for ammonia synthesis reaction. In all these calculations 
but those due to Dyson and Simon [10] either pseudo first 
order kinetics was used or the bulk flow terms in the 
equations for transport within the catalyst pellet were 
asaittod, Dyson and Simon [10] have formifLated the problem by 
considering all the aspects. But unable to solve within 
practicable computation time limits, they also suggest the 
use of an empirical correlation to obtain effectiveness factor. 
The development of transport equations, to follow, are similar 
to those done by Dyson and Simon [10]. 

Catalyst particles are assumed to be spherical [28] 
and isotheimal [30], Knudsen diffusion is neglected [4] and 
diffusion coefficients of each component are assumed to be 
independent of position within a particle. 

The reaction considered for the development of diffusion 
equation is 

2 ^2 2 ^2 


A mole balance for component i gives 



(4.7) 


^2 dr 


■NH, 


oc. 


l-e 


where r is the radial coordinate of spherical particle, 
is the molar fluz of the ith component in the r-direction, 
r^TTi IS the rate of formation of ammonia given hy equation (4.l) 


e IS the void fraction of the packed bed and reprtaeiite the 
stoichiometric coefficient of the ith component in tke 
reaction scheme given by equation (4.6). Following the ordinary 
convention, a^, is taken to be positive if the component is 
a product, negative if the component is a reactant, and zero 
if the component is an inert substance. Nitrogen, hydrogen, 
ammonia, methane and argon have been designated as components 
1 , 2 , 3,4 and 5 respectively. 

Boundary conditions for equation (4*7) are 


f 


dx. 

= 0 at r =0 or, = 0 at r =0 and 
X. = X at r=E’ ( 4 . 8 ) 

where R* is radius of spherical particle, x^ is the mole 

fraction of component i at any point of the catalyst particle 

and X. is that at the surface. It is taken to be same as that 
ig 

in gas. 

The molar fluxes of any two components i and 3 at steady 
state are related as follows; 


a 


i 


N •=» 

] vi 


N. 

1 


(4.9) 
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This implies that molar flux of any inert component 
is equal to zero. The molar flux, of any active component 
i can be expressed in terms of its concentration gradient 
and the molar fluxes of other active components as 


N. = ~C D + X. H. 

T 1 fi . T / I 1 


1 r._l J 
3=1 


(4.10) 


where C is the total concentration of reacting gas mixture 
Kgmole/m^ and D is the effective diffusion coefficient 
of component i. Substitution of equation (4.9) into equation 
(4.10) and utilization of the relationship 


= -1 


(4.11) 


gives 


le dr 


IT = 
1 


(4.12) 


1 + 

“i 


Substituting for IT and IT from equation (4.12) in 

X J 

equation (4.9) yields 


dx. "^i 

-0^3 e-a? 

5 ; 3 = 0^3 [ F 

3 1 


(4.13) 


Integrating this equation subject to boundary conditions 
(4,8), the following equation is obtained 



(4.14) 


4 


a. + X D. /D. 

X = -a +(a + X ) ( ^ 

0 D ' D Dg + 2c^g 

As a consequence of this relationship, the solution 
of only one mole balance (equation 4.7) need be considered. 
Substituting for from equation (4.12) into equation 
(4.7) yields 

dx. 


d^x. 


dr^ 


i. 1 / 1 -s 2 , 


dx 


r dr 


NH. 


a.+x. 

( -i—i) 

/ CD^ J 1-e 
le 


(4.15) 


X* 

Normalising equation (4-. 15) by defining 2 = g*, gives 
d^x. 


dZ 


, dx ^ dx 

1 ( In 2,2 _i 

2 " (a^+x^) dZ'' ^ Z dZ 


NH, 


CD. 

le 1— e 


(4.16) 


and the boundary conditions in terms of Z are 
dx. 

= 0 at Z=Oj X. = X. at z=l (4.17) 

OL 2 1 


This IS a two point boundary value problem and may be 

solved by either choosing tne unknown condition at the centre 

(i.e. x^ at Z=0) and integrating equation (4.16) towards the 

surface (Z=l) for choosing the unknown condition at the 
dx. 

surface (r^ at Z=l) and integrating inward. To solve the 




problem, in both cases, the value of the chosen quantity is' 
varied and integration repeated until the botindary condition 
at the other point is satisfied. It is convenient to write 
equation (4.16) for the product component (EE^), thereby 

- J 

ensuring that the denominator in the second term of eq^mtion 

(4.16) cannot go to zero. If the integration is carried out 

from centre towards the surface, it is necessary to rewrite 

equation (4.16) because of the third term. This type of 

problem has been treated by Weisz and Hicks [38], Hero fhe inward 

integration procedure is used, it is reported [10] that to 

satisfy the boundary conditions at the centre with tolerable 

accuracy it is necessary to determine the flux at the surface 

(i.e. dx^/dZ) to within 1 part in 10^. However, it may not 

be necessary to satisfy the boundary condition at the centre, 

if it IS recognized that the direction of rapid deviation of 

dx^ ... 

the calculated value of from the true one i.e. direction 

dx_ 

of shooting of value as integration proceeds, depends on 

dx_ 

whether the assumed value of Z=l) is smaller or greater 

than the true value. As shown in Figure 4.1, if the assumed 
dx- 

value of is smaller (A) integration would lea'd to shooting 

of its value in one direction as the centre is approached 
whereas it would have the opposite direction if the assuned 
value (B) is greater than the true one. Therefore, by observing 
the direction of shooting, it is possible to obtain a range (AB) 
within which the true value lies. This range could be reduced 




Fig 4 1 - Calculated integratibn path of concentration 
gradient of ammonia in the catalyst pellet 


to desired abcuracy (A*B') By repeated integration, After 

dz, 

aobiieving the desired accuracy (range of value reduced 

to’A'B’) no further integration is called for although the 

"boundary condition at the centre has not been satisfied. 

In this work range A‘B’ representing less than0,5 per cent 

dx, 

of the true value of been used. Mean of the two 

extreme values (A‘ and B’) is taken to represent the true 
dx^ 

values of , 

The surface flux may now be substituted along with 
other terms in equation (4.5) to get the effectiveness 
factor, 


dx, 
dZ '' 


B = 


■R 


! 2 




(4.18) 


The total concentration c is obtained from 
n 


0 = 


S 

1=1 


RT 


(4.19) 


where f^ is the fugacity of component i. The radius, R*, for 
spheres equivalent to industrial size particles may be 
calculated from 

A 


R’ = 


2^ 


(4. 20) 


where TV and ^ are equivalent diameter and shape factor 
respectively. Effective diffusion coefficient is calculated 



( 4 . 21 ) 


m 

from the relation given hy V/heeler [39] 


he = ? ® h 


whore ©■ is the intra -particle porosity and is the hulk 
diffusion coefficient of component i. This coefficient at 
0°G and 1 atmosphere is calculated from the following 
relation 




1 n 


1-z. 




z. 

1 


( 4 . 22 ) 
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where is tne diffusion coefficient of component j in 

component i. The diffusion coefficients calculated from 
equation (4.22) are then corrected for the temperature and 
pressure at the surface of the catalyst pellet by 


^i 



T ^1 .3 

rfy 



( 4 . 25 ) 


where P is in atmospheres. 

4,5 PESCEIPTION OP THE MATHMATIGAL MODEL 

Since the reaction is ezothermic all reactor designs 
have arrangements for removing the heat generated in the 
catalyst bed by the progress of the reaction. The reactor 
also serves as its own heat exchanger to heat the incoming 
synthesis gas. Apart from flow pattern (axial or radial) 
the designs differ on heat removal strategy. In autothermic 
designs, heat is removed throughout the bed volume whereas 



other designs have adiahatic beds with cooling in between 
them either by external exchangers or by quench gas mixing. 
Figure 4.2 shows schematic diagrams of a reactor having 
three adiabatic beds with external cooling of the reaction 
mixture and an autothermal reactor. Both adiabatic as well 
as non-adiabatic beds, in an axial flow reactor, have been 
considered in the present study. 


Adiabatic Catalyst Bed ; 

The temperature and composition throughout a general 
cross-section of the bed is assumed uniform. Axial diffusion 
of mass and heat is neglected. Pressure drop is very small 
compared to the total pressure therefore an uniform pressure 
equal to the average pressure in the reactor has been 
assumed. On the basis of these assumptions the material and 
heat balance equations which describe the evolution of the 


composition and temperature of the system along an adiabatic 


catalytic layer can be written as 

ci 7 


^ ^ H ^^-^3 

dV " ^ G 


(4.24) 


(4.25) 


where G = mass flow rate, 

/XH = heat of reaction, K. cal/lCg.mole 

Y = catalyst volume, m^ and is the extent of 
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reaction defined "by 


= 


'JO 


+ a m 


// 


J 2 - 


<>) 


(4.26) 


where g ,g = mass fraction of component 3 at any point and 
at the inlet, m = molecular weight of component 3 , 

J 


Uon-Adiabatic Catalyst Bed : 

Assumptions same as those in case of adiabatic bed 
are made. Mass balance equation (4.24) holds for this case 
also. The heat balance is described by the following 
equations: 


dT 

S. 

d7 


dV 


U.A.(T-T^) 



Pg 


G 


1 


G. a 


-AE ^ Op 

a mm ■mi ■■ 1 ■■ 1— iw ^ 


a 


G. 


(4.27) 


(4.28) 


where T^ is the temperature of the feed gas in imbedded 
cooling tube, G is the heat transfer coefficient, A is the 
exchange area per unit volume of catalyst and Gpg is the 
specific heat of the feed gas. The equilibrium constant 
K is obtained from the Gillispie and Beattie equation [15]. 
The fugacity coefficients are calculated by means of the 
Cooper expression [9] for nitrogen and hydrogen and by the 
Mewton expression [27] for ammonia. 



The composition and temperature of the reacting system 

A 

at the inlet of any layer being known, the material balance 
equation (4.24) is solved to obtain the extent of reaction 
which in turn gives the mole fraction of individual components 
(equation 4.26). For adiabatic bed the solution of energy 
balance equation (4.25) gives the outlet temperature whereas 
for non-adiabatic bed (cooling tubes imbedded in catalyst bed) 
equations (4.27) and (4.28) are to be solved simultaneously for 
the same purpose, Runge Kutta method of fourth order [8] has 
been used for the simultaneous solution of mass and energy 
balance equations. The size of the integration step has been 
varied according to the volume of catalyst used. 

Calculation Program ; 

On the basis of the mathematical model described above, 
a calculation program in Forxran IV has been prepared for 
IBM 7044 computer. The calculation time for checking the 
performance of a triple-bed reactor with interstage cooling 
is 5 minutes and 3 seconds, when the step size for each bed is 
l/lOOth of the bed volume. 

4.4 RESULTS ARP DISGUSSIOIT 

The calculation results for six different cases having 
adiabatic as well as non-adiabatic beds (autothermal reactors), 
working under different conditions of pressure and composition 
of the reacting system are presented along with the plant data 



in Tables A,1 to 4.4 and Pigures 4.3 to 4.5. 

Por the purpose of discussion the reactors having 
adiabatic bods and those with non-adiabatic beds would 
bo consid>^rod separately. The adiabatic reactors have 
Montecatiri Catalyst w-heroas in autothermal reactors (non- 
adiabatic beds) Haldor Top so e catalyst is in use. 

Re actors with Adiabatic Catalyst Beds : Each reactor under 
consideration has three adiabatic catalyst beds. The results 
for these are presented in Tables 4.1 to 4.3. Pigures 4.3 
and 4.4 show the teiaperature and ammonia per cent profile for 
Case I (Table 4.1). The maxiEiam difference in measured and 
calculated ammonia concentration exists for case I at the 
outlet from first catalyst bed (Table 4.1). The measured 
value IS 4 per cent less than the calculated one. This difference 
IS well within the experimental error limits. Maximum 
difference in temperature is recorded in Case III for exit 
from bed III (Table 4.3). The difference, being only 1.5 
per cent, could be considered insignificant. 

Autothermal Reactors ; Results for these reactors are tabulated 
in Table 4.4. The maximum difference in measured and calculated 
composition exists at the outlet of Case IT (Table 4.4). This 
difference of 3.2 per cent is well within the allowable error 
limits. Moreover, the measured and calculated temperature 
profiles for this case (Pig. 4.5) almost overlap, showing 
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Bed length (fractional) 

Fig. 4- 5 -Temperature of synthesis gas in cooling tubes and catalyst 
bed of an autothermal reactor. 
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complete agreemeiit between the model values and the actual 
oiies, A complete match of temperature profiles woiild mean, 
from theimodynamic considorations, a match in composition 
a,lso, indicating thereby that the difference obse’ived in 
composition may bo because of measurement errors. Also shown 
in Figure 4^5 is the calculatoc' tempera cure profile of the 
feed ga.s in the tubes embedded in the ca,talyst bed. Maximum 
temperature difference of 4 per cent is observed for Case Y, 

For this case the measured and calculated temperature profiles 
show a difference of 2 to 3 per cent over most pert of the 
reactor. In view of vary good agreement for other cases 
this small difference could be considered insignif icanx. 
Moreover, since the comparison is made over the en uire profile 
and not at selected points, the reliability is high inspite 
of small deviations. 

The results very clearly show that in general, the 
calculated values are in very good agreement with the plant 
data,. It IS to be noted that the devia cions, in addition 
to being small, arc randomly distributed with respect to 
actual plant observations both for temperature and composition. 
This confirms the reliability of the rate equations used as 
well as the validitj- of assumptions ma,de under all plant 
conditions. This also shows that the effect of diffusion of 
reactants and products inside the catalyst pores has been 
correctly accounted for. 



MBNCLATUEB 


Constant 

Heat exchange area, 

Total concentration of reacting gas mixture, hg.mole/m^ 
Specific heat of reacting gas mixture, Kcal/kg,®K, 

Specific heat of feed gas, Kcal/kg 

Diffusion coefficient of component i at 0°G and 1 atm. 
Effective diffusion coefficient of component i 
Activation energy for ammonia decomposition, Kcal/lfg,mole 
Eugacity of component i 

“ Eugacity of nitrogen, hydrogen and ammonia respectively 

■ 

Mass fraction of component 3 

Mass fraction of component 3 at inlet 

Mass flow- rate, Kg/hr 

Heat of reaction, Kcal/kg.mole 

Frequency factor in Arrhenius equation for 

Velocity constant of reverse reaction 

Equilibrium constant of the reactionj 0,5 + 1.5 H^M^ 

Molecular weight of component 3 
Molar flux of component i 
Ere s sure, atm. 

Radial coordinate of spherical catalyst particle 
Diffusion coofficiont of componont i 
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NH, 


R* 


T 

T 

U 

Y 

X 


g 


^ig 


Rate of ammonia formation, Kg.mole ’Ml^/hr/m^ of catalyst 
Universal gas constant 
Radius of spherical particle 
Temperature, °K 

Temperature of gas in cooling tubes, 

Heat transfer coefficient, Kcal/hr/m^ 

Catalyst volume, m^ 

Mole fraction of component i 

Mole fraction of component i in gas (hulk phase) 


Greek 

“i 

h 

e 

A 

f 

& 


Symbols : 

Stoichiometric coefficient of component 

Effectiveness factor 

Extent of reaction 

Void fraction of packed bed 

Equivalent diameter 

Shape factor 

Intraparticle porosity 


1 
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CHAPTER 5 


i 

PROCESS SimiATIOI OP AMvIOHIA PLAHT 

There is no recycle of process gas between the point 
of inlet to the PR and the exit from methanator. The outlet 
from one reactor makes the feed to subsequent one. Air 
and steam added at SR and HT inlets respectively are accounted 
for by mass balance. Temperature change associated with 
heat removal in between any two reactors is not considered, 
instead the measured value of temperature after cooling is 
used in simulation calculations. The ammonia synthesis 
reactor involves recycle of a major part of the product gas, 
to feed, after separation of ammonia from it. This recycle 
gas mixed with make up synthesis gas i*e. methanator exit, 

IS the synthesis reactor feed. The process simulation includes 
a scheme for calculating flow and composition of feed to 
ammonia synthesis reactor. The model for this reactor 
calculates the rate of ammonia production. 

In the PR model, dimensions of tube, burner and 
furnace, the number of tubes and burners in use and the 
amount, shape and size of catalyst are required to specify 
the primary reformer. The data on H^/O ratio in hydrocarbon, 
steam, process and fuel hydrocarbon flow rates, inlet 
pressure, inlet temperature, calorific value of fuel, excess 
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air used by burners and the flue gas temperature are required 
by the PR model to calculate the flow,, composition,’ tempera- 
ture and pressure throughout the reformer tube length and at 
the exit. In addition, the total heat absorbed in the 
furnace, the inside and outside tube wall temperatures along 
the entire length are also calculated. 

The outlet from PR is mixed with appropriate amount 
of air at SR inlet, to provide required nitrogen for ammonia 
synthesis and heat for steam reformation of remaining methane. 
The SR simulation model calculates the amount of air required, 
the composition and temperature after its mixing with PR 
exit and the subsequent changes in composition, temperature 
and pressure along the entire length of SR catalyst bed. Only 
data required by the model are dimensions of catalyst filled 
space and amount and specific surface of catalyst pellets. 

The SR exit stream is cooled and in some cases small 
amount of steam is also added before it is fed to the H. T, 

The inlet temperature needs to be specified whereas composition 
of the process gas after mixing with known amount of steam 
IS calculated. Yolume and bulk density of the catalyst 
being known the HT model calculates composition and temperature 
throughout the catalyst volume. 

The HT outlet is the feed to the IT with removal of 
heat in between. The flow and composition remains unaffected 
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but the feed temperature needs to be specified in addition 
to the volume and bulk density of catalyst. The IT model 
calculates composition and temperature throughout the catalyst 
bed. 

The process simizlation assumes complete removal of 
steam from IT exit and complete absorption of CO 2 in the 
CO 2 3-bsorber. It further assumes complete conversion of GO 
in the process gas to GH^ in the mothanator. The exit from 
the methanator 1 , e. the make up synthesis gas therefore 
contains only GH^ and Ar. 

The recycle scheme of ammonia synthesis reactor is 
shown in Figure 5,1. GH^ and Ar which do not take part in 
the reaction are continuously added to the ammonia synthesis 
loop [ABGD] by the make up gas* For continuous operation 
of the process it is, therefore, nocessa3ry that the amount 
of in'orirs being added to the synthesis loop is continuously 
taken out of it. Operating at higher concentration of inerts 
means lower production and higher circulation of synthesis 
gas through the reactor whereas low concentration of inerts 
in recycle gas would require high purge loss. An optimum 
concentration of inerts is therefore maintained. The flow 
and composition of purge, recycle and feed is obtained from 
the following component mass balance equations. 





Mass Bal ance Around Svntliesis Loop (Figure 5,l\ i 


^TSTH 

(5.1) 



M-Xg^ = E-(yg^+ 1.5 yuji^) + 1.5 Prod 

(5.2) 

’^•’=0H^= 

(5,3) 

“■^Ar = E'^Ar 

(5.4) 

Mass Balance Around Point A (Pigure 5,1) • 


M J- R = F 

( 5 . 5 ) 

E-yu^ = P-ffr^ 

( 5.6 ) 

S-yg^ = F-fg^ 

(5.7) 


(5.8) 

M-3c_ 4,. - R-yjir = P-^Ar 

(5.9) 

E-ysHj = E-fHHj 

(5,10) 


Mass Balance Around Poin t B. (Figure 5.l)? 
B = P R + P^od 

B-Zj 5^= (B-i-S)-yu^ 

E-Z^ = (P+R) -y^ 

^2 2 

(r+R)-yoH_^ 

E-Zj^r = (P+E) -y^j. 

= fE+E)-yijH, + 


(5.11) 

( 5 . 12 ) 

(5.15) 

(5.14) 

(5.15) 

(5.16) 


+ Prod 



obtained from equations (5.11) and (5.1^), 

Prod (l“Zjjjj ) 


R 


'Mi, 


^Mi. 


— P 


(5.20) 


Diole fraction of ammonia at the synthesis reactor 

3 

exit is known. Once flow and composition of purge and 

recycle are known, equations (5.5) tnrough (5.10) can be 

solved to obtain the feed flow and composition. However, 

product composition and hence objective of 

3 

synthesis reactor calculation and cannot be known priori. 

Therefore, value is assumed and equations (5.20), (5.5) 

through (5.10) are solved to obtain feed flow and composition. 

The synthesis reactor model calculates the exit concentration, 

' which is compared with its assumed value. In case of 

3 

significant difference a new value for , based on calculated 

3 

value, is chosen and equations (5.20), (5.5) through (5.10) 
solved again to obtain feed conditions. The process is 
repeated till the assumed value agrees with tne calculated 
one within desired limit of accuracy. 

The rate equations used by the synthesis reactor 
model in the present simulation are those for Montedison 
catalyst. In case of a different catalyst it would need to be 
suitably changed. Pressure, volume of catalyst in different 
beds and the inlet temperature for each of them being specified. 



the synthesis reactor model calculates flow, temperature 
and composition of the process gas at each point in the 
reactor and at the exit. The agreement between assumed and 
calculated ammonia concentration in the product gas yields 
the rate of ammonia production. Prom the ammonia cx)ncen- 
tration in the synthesis reacted product gas the rate of ammonia 
production is obtained. The calculation scheme for the 
process simulation is described in detail in Appendix which 
also includes listing of the program. 

Results and Discussion ; 

The calculation results alongwith the measured values 
for two plants at two different levels of operations are 
presented in Tables 5.1 and 5.2. The" 'a~gre^ent- between the 
measured and calculated values of flow, composition, tempera- 
ture and pressure at various point throughout the plant, in 
all the four cases reported, is goncrally vory good* It 
is to be noted that information on process variables required 
by 'the-piant simulation, in addition to process unit specifi- 
cations and inlet conditions, include inlet temperatures to HT 
and LT reactors and to the different beds in ammonia synthesis 
reactor, pressure and the allowable level of inerts in ammonia 
synthesis loop and the per cent of ammonia remaining in the 
gas after separation of the product. The information called 
for by the process simulation is rather very limited whereas 
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TABLE 5.1 

MEAS D'BED AI^TD OALeULAT ED DATA BOB ITAPHTHA BiiSE D 

Aiii' IOBIA PLABT 


Miimber of PR tubes 

Dia of PR tubes 

Heated length of PR tubes 

SR catalyst 

HT Catalyst 

LT ca-talyst 

Synthesis catalyst 

P rocess Inlet PR 

Process naphtha flow, 
kg /hr 

Steam flow, kg /hr 
Pressure, Kg/cm^ 
Temperature, 

SR inlet air temperature 

HT inlet temperature, °0 

LT inlet t emiD3rature,°C 

Synthesis reactor inlet 
t emp era tur e , ° 0 

1st bed 

2nd bed 

5rd bed 

Pressure at synthesis 
reactor inlet, kg/cm^ 


160 

115.3 nm I.B./UT.T 
11,36 meters 
25 tons/l9.4 
75 tons/58.2 
73 tons/56 
57.2 tons/20.9 


Case I 

Case II 

14100 

11900 

67200 

56000 

24.S 

24.3 

443 

440 

C 231 

231 

377 

364 

202 

200 

395 

389 

412 

428 

406 

404 

228 

214 


mm O.D, 
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the measured and calculated values, for different plants, 
compare very well at any point along the flow of process ^s. 
The close agreement in all cases studied shows the versatility 
of the process simulation presented here. 
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MMCLATURE 

E Flow of ^nthesis reactor exit, kg.mole/hr 

F Flow of feed to ^ntiiesis reactor, kg.mole/hr 

f^ Mole fraction of component i in feed 

P Flow of purge gas, kg.mole/hr 

Prod Rate of a3iimonia production, kg.mole/hr 

R Flow of recycle, kg.mole/hr 

x^ Mole fraction of component i in make up synthesis gas 

y^ Mole fraction of component i in purge and recycle gas 

Mole fraction of component i in synthesis reactor 
exit 

Chemical formula is used as subscript to represent 
the respective component molecules. 



GBABTEh. 6 


CQJCLUSIOJG AJJ Pu:.C0i^-I5ITI}.x;iQ::S 

Mathematical models for simulation of all the important 
process units of ammonia plant, namely, primary and secondary 
reformers, high temperature and low temperature shift 
reactors and ammonia synthesis reactor, have been developed. 
These have been subseq_uently used for process simulation of 
an ammonia synthesis plant. 

In primary reformer model, rate equations for steam 
methane reactions and heat transfer equations have been 
successfully developed and used. Transfer of heat from 
flames which accounts for 25-35 per cent of the total heat 
transfer has separately been considered. 

Secondary refoimer model calculates the amount of air 
to be mixed, with PK exit, at the inlet to SR, the resulting 
temperature of process gas and the subsequent; change in 
composition and temperature over the catalyst bed. Rate of 
reformation reaction in the SR is obtained by considering 
the reaction to be limited by the rate of diffusion of 
reactants from bulh to the catalyst pellet surface, 

A rate equation for water-gas shift reaction over 
chromia promoted iron— oxide (HT catalystj has been success- 
fully developed and used in simulation of several industrial 
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reactors. Rate of shift reaction over IT catalyst (CuO-ZnO) 
has also been aeveloped and used. 

A suitable honiogenoous rate expression for anDionia 
synthesis reaction has been selected, A iiethod of calculating 
effectiveness factor for industrial catalysts, by partially 
solving the transport equations has been developed. Subse- 
quently, these have been used in a laatheinatical model developed 
for design and simulation calculations of axial flow ammonia 
synthesis reactors. Reactors having adiabatic as well as 
those With non-adiabatic beds have been considered. 

The process simulation of auuonia plant has been obtained 
using these models and considering other process units to 
be ideal. Recycle of gas at the ^nthesis reactor is also 
considered. The calculation results and the measured values 
of process variables at different locations in a plant are 
in close agreement. 

The models developed in this work have been used only 
to check the performance of various units but it can be claimod 
that their use for design work would yield very good results. 

Recommendations ; 

Efforts are required to develope methods to calculate 
transfer of heat to RR tubes whrtk the flue gas temperature 


IS not uniform. 
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E,a‘ce equations developed hero for wator gas shift 
reaction over HT and IT catalysts and auuonia synthesis 
reaction over iron catalyst are valid for the catalysts being 
used. Moro extensive work is required to develop truely 
general rate equations which will be valid for any catalyst. 

Attempts should be nade to develop such a siciulation 
ol aniionia plant which would account for the perf oruanco of 
all the units and also consider the recovery of heat fron 
SR exit, PR exit, flue and synthesis gas and power generated 
for coupression of synthesis gas. 
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APPEiroiX 


PROC E SS SIMLATIOE PROGRAM POR AlfflOSlA PIAR T 


The process simulation program consists of a mam 
program and fifteen subprograms, namely PIREBX, BQMC, ¥ALL, 
COMP, SRIITC, SHIRT AIMORIA, STOP, REACR, REAEH, EPP, EPPFH, 
HTKCH, CP and RIJICE. 

in Program ; 

All relevant data about the plant is obtained by this 
program through READ stabements. These data are printed in 
the same format to enable a check on the possible errors. 

The required data are transferred to the subprograms whenever, 
they are called in bjr the main program, through the arguments 
and COMMON statements and in tan these are transferred as 
required to the subprograms that these subprograms call and 
so on. 

After reading and printing data the main program calls 
PIREBX to calculate emissivities of flue gas and flemes and 
the total amount of heat absorbed in reforming furnace. These 
values along with the flow and composition of flue gas are 
printed by xhe subprogram. Composition of PR inlet stream 
is calculated in main program by assuming conversion of feed 
to methane and carbon dioxide only. Subsequenbly SgiC is 
called to calculate the real composition of feed containing 
CH^, CO, CO^ and Plow, composition, temperature and 





pressure along "the flow over the catalyst in the PE tube 
are caleulated by the integration subprogram RUNGS whach 
calls in WALL to ealeulate inside and outside tube wall 
temperatures, ESACE to calculate rates of steam-methane and 
water-gas shift reactions, HIEXN to calculate heat of the 
two reactions, CP to calculate specific heat of reaction 
mixture, COMP to calculate the composition of the process 
gas and HTCP to calculate inside tube heat transfer coefficient. 
Temperature, component flow rates and pressure of the process 
gas, outside and inside tube wall temperatures and amount of 
heat transferred to a differential section of the reformer 
tube length are calculated and printed at each integration 
step. 

SRINTC IS called by the mam program to calculate the 
amount of air req.uired and the temperature and composition of 
the mixed feed. EUNGh is again called to integrate the 
changes in flow, composition, temperature and pressure over the 
SE catalyst bed. Rate of steam-methane reaction is calculated 
from diffusivity considerations which involves equilibrium 
calculations to be done by EQi^IC. EEICE is called to calculate 
rate of shift reaction. Heat of reactions, specific heat of 
process gas and composition of the reaction product at any 
point along the path of integration are evaluated by HTRXN, 

CP and COMP respectively. Volimie of the catalyst and correspond^ 
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"beniperature , pressure and flow rates of each, component of 
process gas are pf-inted for the entire volume of the 
cat alyst . 

Any steam added to the process gas at the inlet is 
accounted for to obtain the inlet process gas composition 
and this infomation along with inlet temperatures, pressures 
and catalyst volumes for HT and If reactors are transferred 
to SHIFT which is called in to calculate the performance of 
HT and IT reactors, SHIFT calls in ROTGE to perform inte- 
gration Ox mass and heat balance equations whereas HHACR, 
HTRIH and CP are called to calculate reaction rate, heat of 
reaction and specific heat respectively. Volume of the 
catalyst alongwith the calculated temperature and composition 
of the gas throughout the entire volume of HT and IT reactors 
are printed. 

The CO2 and H2O content are taken as zero. The compo- 
sition change associated with the removal of CO2 and H2O and 
complete conversion of CO to CH^ is calculated. This gives 
the make up synthesis gas composition. Flow and composition 
of recycle gas is calculated by assuming the ammonia concen- 
tration in synthesis reactor outlet gas. The composition of 
feed gas i.e. the mixture of make up and recycle synthesis 
gas IS calculated and this along with the volume of the 
catalyst in each bed and the corresponding inlet temperature 
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is traiisf jrrod to ildOFlA. This subprogram calls EPFEE to 
calculate the of f octivc!n''Ss factor of. the synthesis cataljrst 
at each step oi intogration which is don: by RUFO-S, EFENH 
cannot bo called by RblTGE since this itself calls RUFGE, 

Therefor: for th"' four subst^ps tahon by 3Uh&3 to complete 
one intogration step, the offoctivcnss faexor is taken 
constant. EUHGE calls RBiUTH, KTEXh and CP to calculate the 
rare of ammoma synthesis reaction, heat of reaction and 
specific heat of synth^^sis gas respoctivc'ly. Effectiveness 
factor, voluiae of thL catalyst, temperature and composition of 
gas are printed at each integration stop. If tho calculated 
ammonia concentration of the produce gas does not agree with 
th'‘ assumed value, calcula/cD-oxi is rop’''atcd mth a new value of 
ammonia concentration in product gas. In case of agreement 
between tho assumed and calculabed value of ammonia concentration 
at the sjrnthosis reactor ourlot, tb^ rate of ammonia production 
IS evaluated and printed. e 

S ubprogram ; 

pip-FjP.X ; Data required by this subprogram arc - Hydrocarbon 
fuel used, H 2 /C leatio, flow rate- and calorific value of tho 
fuel, beam lengths of flames and flue gas, ambient t omp ^.ra bui e , 
excess air used in burning tho fuel and temperature of flue 
gas. It calculates the adiabatic flame temporature by 
equating sum of tho calorific heat in furl and the sensible 
heats in flue and flames to the sensible heat in flames. 







Amount of lieat absorbed in reformer is the difference between 
sensible heat in flames and that in flue gas leaving the t 

reformer. Adiabatic flame temperature, emissivities of flue 
and flames, and the total heat absorbed m primary reformer 
are printed, 

BQMQ: lata required are component molar flow rates, pressure 
and temperature. Equilibrium composition at the given 
pressure and temperature is calculated by ropoatodly evaluating 
mole fractions of and CO from equilibrium of steam-methane 
and water-gas shift reactions respectively. Calculated 
equilibrium composition is transferred to the point which - 
called the subprogram and no values are printed. 

WALL : In addition to the reformer specifications, data 
required are the temperatures and emissivities of flue gas 
and flames, mass flow rate, pressure, temperature and 
composition of process gas and specific heat of reaction 
mixture. Inside and outside wall temperatures, heat transferred 
from flue gas and flames to the differential section under 
consideration are obtained oy equating tho amount of heat 
transferred from flames and fine gas to that from inner .wall 
of the tube to the process gas. Wo values are printed. 

COMP : Stoichiometric coefficients, molecular weights and mass 
fractions of all component molecules, total mass flow rate 
and extents of reactions are the data required. Composition 



of the process gas for given extent of reactions is calculated 
hy the subprogram. 


SRINTO : Molecular weight, temperature and molar flov^ rates 
of different components in IR exit and temperature of air is 
to be provided for tnis subprogram to calculate amount of 
air required to be mixed t* process gas and resulting tempera- 
ture, composition and flow. 

SHIJ'T % Stoichiometric coefficients for shift reaction and 
molar flow rates of ditleient comjjoaenc;s at tne inlet to hi, 
pressure, volume of catalysts in h'T and IT and the inlet 
temperatures for the two reacuors are tho information called 
for. Mass and energy oalance equatxons over a drfferentral 
section of reactor are intet,rated over the entire volume, 
from ixilet to exit, with RuRGn. The reaction rates over the 
HT and the LT catalysts are colcolated by RBACR. Effective- 
ness factor is obtained from ERR, heat of reaction from HTRZR 
and specific heat from OR. At each integration step, volume 
of catalyst, temperature and composition of the process gas 
are printed. 

Data required include stoichiometric coefficients 
for synthesis reaction, molecular weights, molar flow rates 
of all components and volume of catalysb ano. inlet temperatures 
for all the catalyst beds. ElFHH is calleo. to calculate 
effectiveness factor lor the catalyst and subsequently RURGE 



xs called to integrate the mass and heat balance equations 
over a differential section which needs KBiJH, HTRItJ and 
CP to calculate rate of synthesis reaction, heat of reaction 
and specific heat of the process gas respectively. Volume 
of the catalyst, temperature and composition of process gas 
and effectiveness factor are printed at each integration 
st ep. 

HT Cr' i Diameter of the catalyst pellet, specific heat, temperature 
and process gas mass flow rate per unit area aro supplied to 

this subprogram for calculating reformer tube internal heat 
transfer coefficient. 

REAGR; Composition, temperature and pressure of the process 
gas along with the catalyst under consideration being specified 
this subprogram calculates the rate of methane-steam reaction 
over reformation catalyst and that of shift reaction over 
reformation, HT and IT catalysts. 

CP; This subprogram calculates the specific heat of process 
gas throughout the ammonia plant which has components from 
among the eight considered by this subprogram. Composition, 
temperature and pressure being specified, this subprogram 
calculates the average molar specific heat of process gas. 

BJBAM; This subprogram calculates rate of ammonia synthesis 
reaction from component gas fugacities and temperature of 


process gas. 
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EFg ? Effectiveness factors for HI and IT catalysts are 
obtained from this subprogram given the forward rate constant, 
temperature of reaction and the reactor under consideration. 

EPF HE ; Pugacity coefficients, bulk dif fusivities and mole 
fractions of component molecules, stoichiometric coefficients 
for synthesis reaction, temperature and pressure of the gas 
are the data required. The subprogram calculates effective- 
ness factor for the ammonia synthesis catalyst particles. 

HTRZH ; Pressure and temperature of process gas and the 
reaction under consideration being specified, this subprogram 
calculates the heat of tne specified reaction. Heat of 
reformaxion, shift and synthesis reactions, one at a time, 
can be obtaaned from this subprogram. 

RUNG-E ; 1 integrations involved in the program are performed 

by tnis subprogram. This employs the fourth order Runge-Kutta 
method witn Eutta' s coefficients to integrate a system of H 
simultaneous first order ordinary differential equations 

F(J) = DT(J)/M, (J=l,2, ,H) 

accross one step of length H in the independent variable Z, 
subject to initial conditions Y( J) , (J = 1,2,...,H), N, H, P 
and initial values of Z and Y(J), ( J=l, 2 , . . . ,H) are required 
to perform the integration. 
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